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QQ ' Abstract. By combining our UVES-VLT spectra of a sample of four damped Lya systems (DLAs) toward the 

quasars QOlOO-l-13, Q1331-I-17, Q2231-00 and Q2343-I-12 with the existing HIRES-Keck spectra, we covered the 
total optical spectral range from 3150 to 10 000 A for the four quasars. This large wavelength coverage and the 
^ I high quality of the spectra allowed us to measure the column densities of up to 21 ions, namely of 15 elements — N, 

^ . O, Mg, Al, Si, P, S, CI, Ar, Ti, Cr, Mn, Fe, Ni, Zn. This comprehensive set of ionic column densities and elemental 

abundances severely contrasts with the majority of DLAs for which only a handful of ions and elements is typically 
observed. Such a large amount of information is necessary to constrain the photoionization and dust depletion 
effects, two important steps in order to derive the intrinsic chemical abundance patterns of DLAs. We evaluated 
the photoionization effects with the help of the Al^/Al^""", Fe^/Fe^"'", N^/N^ and Ar/Si,S ratios, and computed 
, dust corrections. Our analysis revealed that the DLA toward Q2343+12 requires important ionization corrections. 

This makes the abundance determinations in this object uncertain. The access to the complete series of relatively 
robust intrinsic elemental abundances in the other three DLAs allowed us to constrain their star formation history, 
their age and their star formation rate by a detailed comparison with a grid of chemical evolution models for 
Q ' spiral and dwarf irregular galaxies. Our results show that the galaxies associated with these three DLAs in the 

^ , redshift interval z^ha = 1.7 — 2.5 are either outer regions of spiral disks (radius > 8 kpc) or dwarf irregular galaxies 

C/3 ■ (showing a bursting or continuous star formation history) with ages varying from some 50 Myr only to > 3.5 Gyr 

^ ' and with moderate star formation rates per unit area of —2.1 < logi/; < —1.5 Mq yr~^ kpc~^. 
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1. Introduction chemical enrichment and ionization histories of gaseous 

structures on a variety of scales, ranging from the inter- 
The study of high redshift galaxies and the access to their galactic medium to high column density absorption sys- 
physical properties can be done using either emission or ^^^^ associated with galaxies. Third, since the detection 
absorption Ime spectroscopy. The absorption line spec- material intercepting a line of sight to a given QSO is 
troscopy is a very powerful technique, and it presents sev- dependent only on the column density of the gas and the 
eral advantages compared to the information which can i^j^inosity of the QSO, this is a unique technique for prob- 
be derived from emission. First, it allows to detect objects -^^ chemical composition and physical conditions in 
up to very high redshift. Indeed, the quasars (QSOs) de- ^j^^ interstellar medium of various types of galaxies over 
tected up to redshifts of 6.2 (" Fan et alj | 200lD can be used ^ j^^g^ ^^^g^ lookback times, detected independently 
as background searchlights for useful probes of the inter- ^j^^j^ distance, luminosity, star formation history and 
veiling Universe up to lookback times of 95% of the age of morphology 
the Universe. Second, through the analysis of QSO absorp- 
tion lines we can study the spatial distribution, motion. While the major^ part of the baryon content of the 

present-day galaxies is concentrated in stars, it must exist 

Send offprint requests to: M. Dessauges-Zavadsky an epoch at which the essential of the mass of galaxies 

* Based on the UVES observations collected during the still lay in the gas. The damped Lya systems (DLAs) ob- 

ESO programmes ID No. 65.O-0296 and 67.A-0022 at the served in QSO spectra are characterized by large column 

VLT/Kueyen telescope, Paranal, Chile densities of neutral hydrogen, > 2 x lO^'' cm~^, and domi- 
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nate the cosmic mass density o f neutral hydrogen gas (e.g. 
Storrie-Lombardi et a l. 1996; IStorrie-L ombardi & Wolfg 
200(11 ■ The general view is that the DLA systems probably 
represent some early stages in the evolution of the galaxies 
we see around us today, perhaps at a time shortly after 
they had condensed out of the intergalactic medium, but 
before they had time to form many stars, so that most of 
their mass still resided in the interstellar medium. These 
objects thus are by far our best laboratory for studying the 
galaxies at high redshift, in their early stages of evolution, 
and for tracking the galactic chemical evolution through 
the cosmic ages. 

However, the reconstruction of the star formation his- 
tories of DLAs from the abundance pattern measurements 
is not s traightfo rward and led to contradictory results 
iILu et al...l99&:.Prochaska fc Wolfelll999HCenturi6n et all 
l2000t iMolaro et alJl2000t IVladild 120*0231) . The access to 
abundance ratios involving two elements formed on dif- 
ferent timescales, in particular the a/Fe ratios, examined 
together with [Fe/H], or any other metallicity tracer such 
as [Zn/H], is crucial, since the star formation history of 
a galaxy is completely d etermined by th e [a/Fe] versus 
[Fe/H] distributions Ce.g. lMatteuccill200l|) . The principal 
difficulty, however, is to disentangle the nucleosynthetic 
contributions from dust depletion effects. Because we are 
studying gas-phase elemental abundances in DLAs, the 
observed abundances may not represent the intrinsic com- 
position of the system if part of the elements is re moved 
from the gas to the solid phase ijSavage fc Sembachl fl996') . 
Several pieces of evid ence show that some dust is indeed 
prese nt in DLAs (e.g. IPei et a,ljll99ll IProcha,ska, fc Wolf j 
l2002t) . Another effect which has to be carefully examined, 
although usually assumed negl igible is the p hotoionization 
effect (e.g. .Viegaa. 1995: Vladilo et aLlbOoH) . since in gas- 
phase studies a fraction of the gas may also be ionized. 

The difficulty to correctly evaluate the photoionization 
effects and the degeneracy between dust depletion and nu- 
cleosynthesis in DLAs is accentuated by the limited num- 
ber of ions and elements typically detected in these galax- 
ies. To determine the "intrinsic" chemical abundance pat- 
terns, free from ionization and dust depletion effects and 
to be then able to study the DLA galaxies individually, we 
need to examine several column density ratios of adjacent 
ions of the same element and the relative abundances of 
as many elements as possible. We thus aimed at obtain- 
ing for a few damped Lya systems the column density 
measurements of many ions of different ionization levels 
and the abundance measurements of the complete series 
of accessible elements in DLAs. Until now the DLA galaxy 
population has been analyzed as a whole and chemical 
evolution models were constructed in order to interpret 
the abundance patterns observed in DLAs as an ensem- 
ble, considering them as an evolutionary sequence (e.g. 



Matteucci et 31.11997": 'Jimenez. 


Bowen fc Matteuccll999l 


Hou et al.ll200lU'Mathhn et al. 


200ll 


, while several pieces 



of evidence — the low redshift deep imaging revealing a 
variety of morphological types belonging to the DLA pop- 
ulation (e.g. .Le Brun et al... 1997.: .Nestor et al...2002) . the 



large scatter in the a over Fe-peak element abundance ra- 
tios at a given metallicity and the large scatter observed in 
the metallicitics — show that the DLAs trace galaxies with 
different evolutionary histories. Some may have formed 
stars on timescales similar to that of the early Milky Way, 
while others apparently did so more slowly or intermit- 
tently, so that the Fe-peak elements could catch up with 
the Qf-elements. The DLA systems thus likely sample a 
wide range of galaxy types, and consequently a variety of 
star formation histories. It is very important to determine 
the star formation history of each of these high redshift 
galaxies individually to better understand the galaxy for- 
mation and evolution. 

We thus attempted a new approach for studying the 
DLA galaxy population, focusing on individual systems. 
By combining our UVES-VLT spectra with the existing 
HIRES-Keck spectra we obtained the column density mea- 
surements of some 21 ions, namely of 15 elements — N, O, 
Mg, Al, Si, P, S, CI, Ar, Ti, Cr, Mn, Fe, Ni, Zn in four 
DLA systems in the redshift interval ^abs = 1.7 — 2.5. The 
complete analysis of the four DLA systems is described 
in Sects. 121 and 121 In Sect. ^ we address the question of 
photoionization effects, and in Sect. Owe analyze the dust 
content and discuss for each system individually the intrin- 
sic chemical abundances. Finally, the relative abundances 
as a function of the metallicity and the redshift are sys- 
tematically compared with a grid of chemical evolution 
models for spiral and dwarf irregular galaxies in Sect.O 

2. Observations and data reduction 

The selected quasars QOlOO-l-13, Q1331-I-17, Q2231-00 
and Q2343-I-12 with four intervening DLAs in the red- 
shift interval Zabs = 1.7 — 2.5 are relatively bright 
with V=16.5— 17.5 and already known from the liter- 
ature. Indeed, the DLAs toward Q0100-fl3, Q1331-hl7 
an d Q2231— 00 were c areful ly a nd accurately analyzed 
bv lProchaska fc Wolfe! l|l999l) and lProchaska et alJ (|20o3) 
thanks to high-resolution spectra obtained with the 
HIRES echelle spectrograph on the Keck I telescope at 
Mauna Kea, in Hawaii, and some results also obtained 
from HIRES-Keck spe ctra on th e DLA toward Q2343-I-12 
were presented by e.g. lLu et alJ l)l998j) . 

To complete the wavelength coverage of the HIRES- 
Keck spectra of these quasars, we used the unique capabil- 
itv of the Ultra violet- Visual Echelle Spectrograph UVES 
((D'Odorico et al. 2000) on the VLT 8.2m Kueyen ESQ 
telescope at Paranal, Chile, and we obtained high resolu- 
tion, high signal-to-noise ratio spectra for each quasar in 
the blue A 3150-4500 A and in the far-red A 6700-10000 
A. The observations were performed in visitor mode in 
September 2000 for one object and in service mode in 
summer 2001 for the three other objects under good see- 
ing conditions (between 0.5" and 1.0"). For each obser- 
vation, slit widths of 1" in the blue and of 0.9" in the 
red were chosen with a CCD binning of 2 x 2 resulting 
in a resolution of FWHM ~ 6.9 km s^^ and 6.4 km s^^, 
respectively. Relevant details of the observations are col- 
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Quasar 


V 


2cm 


Zahs 


Observing date 


Mode 


Wavelength range 


Exposure time 




[mag] 










[nm] 




[s] 


QOlOO+13 


16.6 


2.69 


2.309 


September 2001 


dic2 (390+860) 


330-452/670- 


-1000 


7200 


Q1331+17 


16.7 


2.08 


1.776 


June- July 2001 


did (346+860) 


312-387/670- 


-1000 


13500 


Q2231-00 


17.5 


3.02 


2.066 


September 2000 


dic2 (390+800) 


335-452/610' 


-990 


15700 


Q2343+12 


17.0 


2.55 


2.431 


August 2001 


did (346+860) 


315-387/670- 


-1000 


14400 



lected in Tabled The total exposure times of each quasar 
were split in multiple exposures of 3600 or 4500 s. 

The spectra were reduced using the ESO data reduc- 
tion package MIDAS and the UVES pipeline in an inter- 
active mode available as a MIDAS context. A detailed de- 
script ion of the pipeline can be found in iBallester et alJ 
l|2000t> . To make sure of the best result, we made a sys- 
tematic check of each step of the pipeline reduction. The 
wavelengths of the reduced one-dimensional spectra were 
converted to a vacuum-heliocentric scale, and the indi- 
vidual spectra of each object were co-added using their 
signal-to-noise ratio as weights in order to get the maxi- 
mum signal-to-noise ratios. The final step was the normal- 
ization of the resulting spectra obtained by dividing them 
by a spline function fitted to smoothly connect the regions 
free from absorption features. The continuum in the Lya 
forest was fitted by using small regions deemed to be free 
of absorptions and by interpolating between these regions 
with a spline. An average signal-to-noise ratio per pixel 
of ^ 30, 55 and 45 was achieved in the final spectra at A 
- 3700 A, 7500 A and 9000 A, respectively. 



3. Data analysis and ionic column densities 

By combining our UVES-VLT spec tra with the exist- 
ing HI RES-Keck spectra obtained bv lProchaska fc Wolfel 
l)l999(l we cover the total spectral range from 3150 to 
10000 A for the four observed quasars, and hence have 
access for the first time to up to 21 ions and 15 elements 
for each of their intervening DLAs. 

The column densities of the metal species were derived 
with the Voigt profile fitting technique. This technique 
consists in fitting theoretical Voigt profiles to the observed 
DLA absorption metal lines well described as a complex 
of components, each defined by a redshift z, a Doppler 
parameter b, a column density N and the corresponding 
errors. The fits were performed u sing a minimization 
routine fitlyman in MIDAS (Fontana fc Ball ester 1995). 
We assumed that metal species with similar ionization 
potentials can be fitted using identical component fitting 
parameters, i.e. the same b (which means that macrotur- 
bulent motions dominate over thermal broadening) and 
the same z in the same component, and allowing for vari- 
ations from metal species to metal species in N only. We 
distinguish three categories of metal species with simi- 
lar ionization potentials: the low-ion transitions (i.e. the 
neutral and singly ionized species), the intermediate-ion 



transitions (e.g. Feiii, Aim), and the high-ion transitions 
(e.g. Civ, Siiv). We used relatively strong (but not sat- 
urated) lines to fix the component fitting parameters (5 
and z), and excellent profile fits could then be achieved 
for weak metal lines and for metal lines located in the 
Lya forest where the probability of blending is high by 
allowing only the column density to vary. We had a suffi- 
cient number of relatively strong metal- line profiles to well 
constrain the fitting parameters in the four studied DLAs 
exhibiting multicomponent velocity structures. 

Throughout the analysis we adopted the list of atomic 
data — laboratory wavelengths and oscillator strengths — 
presented in Table [3 The most recent measurements of 
A and /-values of the metal-ions that impact the abun- 
dances of DLA systems and their references are reported 
there. Information on the at omic data of add i tional use- 
ful metal-ions can be found in lProchaska et al and 
on the web site of "The HIRES Damped Lya Abundance 
Database"^ maintained by J. X. Prochaska and coUabo- 
rator s. Finally, we adopt the so lar meteoritic abundances 
from lGrev esse fc Sauvall l|l998l) . 

The measured component per component ionic col- 
umn densities obtained from the fitting model solutions 
of the low- and intermediate-ion transitions are summa- 
rized in Tables ISJ-Q The reported errors are the la errors 
on the fits computed by fitlyman. They possibly under- 
estimate the real error on the measure, since they do not 
take into account the uncertainty on the continuum level 
determination. For the components where the line pro- 
file is saturated, the column densities are listed as lower 
limits. The values reported as upper limits are cases with 
significant line blending with H i clouds of the Lya forest 
or with telluric lines. The fitting solutions of the low- and 
intermediate-ion transitions are shown in Figs. Eland 
IHlfor the four DLA. In these velocity plots, v = corre- 
sponds to an arbitrary component, and all the identified 
components are marked by small vertical bars. The thin 
solid line represent the best-fitting solution. The telluric 
lines have been identified thanks to the spectra of a hot, 
fast rotating star taken in the same nights as the scientific 
exposures. 

The neutral hydrogen column densities were estimated 
from the fits of the Lya damping line profiles. The 6- values 
were fixed at 20 km s^^, and the redshift z were left as 
a free parameter or fixed to the redshift of the strongest 
component of the metal-line profiles according to the sys- 
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Table 2. Atomic data 
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0.07786000 
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tern (see comments in the following sub-Sections). When 
other lines of the Lyman series were accessible in our spec- 
tra, we used them to check the H I column densities derived 
from the Lya lines. Figures [S] and [7| show the results 
of the H I fitting solutions for the four DLA systems. The 
small vertical bar corresponds to the redshift obtained in 



the best-fitting solution and the thin solid line represents 
the best fit. 

We now briefly describe the fitting results for the four 
individual DLA systems. Their Mnii A2576,2594,2606, 
Till Arost > 3000 and Mgii A1239,1240 fits and col- 
umn density measurements have already been discussed 
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Table 3. Component structure of the Zabs = 2.309 DLA 
system toward QOlOO+13 



No Zabs 






Ion 


logAf(o-iogiv) 






km s~^ 


km s~^ 






Low-ion transitions 


1 


2.309027 





5.6(0.1) 


Fell 


14.86(0.01) 










Znll 


12.28(0.01) 










Cm 


13.15(0.01) 










Sii 


14.85(0.05) 










An 


13.99(0.15) 










Pii 


12.88(0.08) 










Ni 


14.79(0.08) 










Mgii 


15.41(0.07) 


2 


2.309161 


+ 12 


5.3(0.1) 


Fell 


14.61(0.01) 










Znll 


12.03(0.02) 










Cm 


12.96(0.01) 










Sii 


14.73(0.08) 










An 


13.81(0.08) 










Pii 


12.55(0.10) 










Ni 


14.65(0.12) 










Mgii 


15.06(0.13) 


3 


2.309443 


-1-37 


14.1(0.2) 


Fell 


14.02(0.03) 


4 


2.309486 


+52 


2.8(0.3) 


Fell 


12.96(0.03) 


Intermediate-ion transitions 


1 


2.308980 


-4 


6.0(1.9) 


Aim 


12.18(0.12) 










Feiii 


< 13.28 










Nil 


< 13.24 


2 


2.309093 


+6 


11.2(1.2) 


Aim 


12.22(0.08) 










Fclll 


< 13.22 










Nil 


< 13.64 


3 


2.309498 


+43 


13.0(1.9) 


Aim 


12.12(0.05) 



Velocity relative to 2 = 2.309027 



in details in iDessauees-Zavadskv et all l)2002al hereafter 
Paper I). Let's just remind that the Mgii A1239,1240 hues 
are detected in the DLA Lya damping hne red wing and 
that they have been fitted after a local renormalization of 
the spectrum around the Mg ii lines with the fit of the Lya 
damping wing profile. The errors on the measured Mg+ 
column densities have been estimated by varying the con- 
tinuum level by 5%. 



3.1. QOlOO+13, z, 



abs 



2.309 



and 



This system was careful ly studied bv IWolfe et al.1 
and subsequent l y by IProchaska 
IProchaska et alJ l|200lh . Thanks to the UVES spectra 
we confirm some of their column density measurements, 
namely A^(Zn+), A^(Cr+) and A^(Fe+) obtained from the 
Fell A2249,2260,2344,2374 lines not present in the HIRES 
spectra, and we obtain new column density measurements 
of Sii, An, Pii, Ni and Mgii (see Fig. E)). We also get 
upper limits on the column densities of two intermediate- 
ions. Nil and Feiii located in the Lya forest, by us- 
ing the profile shapes of the intermediate-ion lines Aim 
A1854,1862 observed in the HIRES spectra to constrain 
the fitting parameters (see Fig.EJ. 
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4150 
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3360 3380 3400 3420 

Observed Wavelength [A] 

Fig. 1. Normahzed UVES spectra of QOlOO+13 showing 
the DLA Lya and Ly f3 profiles with the Voigt profile fits. 
The vertical bar corresponds to the wavelength centroid 
of the component used for the best fit, z — 2.309027. The 
measured Hi column density is logA^(Hi) — 21.37 + 0.08. 



The low-ion absorption line profiles of this system are 
characterized by 4 components. Their redshifts, ^-values 
and column densities are presented in Table |3| Two of the 
4 components, the components 3 and 4, are too weak to 
be accurately identified in the weak metal-lines and in the 
metal-lines located in the Lya forest. Thus, we prefer to 
avoid measuring them in most of the observed transitions. 
Being weak and since we are interested in the relative el- 
emental abundances, this has no impact on the final con- 
clusions on this DLA system. The intermediate-ion lines 
show very similar profiles to the low-ion profiles suggesting 
they are coming from the same absorption regions. Their 
fitting solutions are given in Table |3| 

Figure^shows the fitting solutions of two Lyman lines 
of this DLA, Lya and Ly/3, observed in the UVES spectra. 
The fits were obtained by fixing the 6-value at 20 km s~^ 
and the redshift at z = 2.309027, i.e. at the redshift of the 
strongest metal- line component (the component 1). The 
measured Hi column density, logA^( Hi) = 21.37 + 0.08 , 
is very similar to the one derived bv lPettini erall l|l990tl 
from low-resolution spectra. Giving the high H° column 
density, the red wing of the Lya damping profile extends 
over the Sii A1250,1253,1259 lines. Similarly to the Mgii 
A1240 hne (the Mgii A1239 hne being heavily blended), 
we had first to locally renormalize the spectrum with the 
fit of the Lya damping wing profile before measuring the 
S+ column density. 



3.2. Q1331+17, Zabs = 1-776 

This famous DLA system has been studied by a large 
number of authors . We m entio n here only the work o f 
IProchaska fc Wolfd ((l£9^ and IProchaska eT^ l|200lh 
who obtained the highest quality data and made the most 
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Fig. 2. Velocity plots of the metal line transitions (normalized intensities shown by dots with la error bars) for the 
DLA toward QOlOO+13. The zero velocity is fixed at z = 2.309027. For this and all the following figures with velocity 
plots, the vertical bars mark the positions of the fitted velocity components and the symbols © correspond to the 
telluric lines. The letter K refers to the Keck spectra. 



accurate analysis. By extending their wavelength cover- 
age with our UVES data, we obtained additional col- 
umn density measuremen ts of Sii, Pii, Ni, Mgii, Mnii 
and CI I. We confirm t he IProchaska fc Wolfel l)l999() and 
IProchaska et alJ l)200l[) column density measurements of 
Sill, Fell, Ci, Mgi and Aim (see Fig. But, their 



Af(Ni+) and A^(Cr+) measurements differ from ours by al- 
m ost 0.1 dex. We a lso confirm the revised value of A^(Zn+) 
of IProchaska et aLi (2001) obtained by correcting the con- 
tamination of Znii A2026 by the Mgi A2026 profile. The 
detection of the relatively strong Mgi A2852 line in the 
UVES spectra has allowed us to measure even more reli- 
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Table 4. Component structure of the Zabs = 1.776 DLA system toward Q1331+17 



No 2abs 






Ion 


log7V(criogjv) 


No 


^ahs 


"rcl 




Ion 


logAf(o-iog]v) 






km s^^ 


km s^^ 










km s~^ 


km s~^ 






Low-ion transitions 




1 


1.776336 


-4 


12.0(0.3) 


Fell 


13.90(0.04) 


4 


1.776724 


+38 


13.4(0.8) 


Fen 


14.07(0.02) 










Cm 


9.00(0.25) 










Znll 


11.50(0.06) 










Niii 


12.82(0.13) 










Cm 


12.22(0.06) 










Mnii 


11.84(0.06) 










Nin 


12.87(0.04) 










Mgii 


11.00(0.25) 










Mnii 


11.91(0.01) 


2 


1.776370 





9.8(0.1) 


Fell 


14.16(0.02) 










Sin 


14.60(0.01) 










Znii 


12.35(0.01) 










Sn 


14.15(0.11) 










Cm 


12.67(0.01) 










Pn 


12.41(0.21) 










Niii 


12.89(0.10) 










Ni 


< 14.30 










Mnii 


12.04(0.03) 










Mgn 


14.99(0.14) 










Sin 


15.06(0.01) 


5 


1.776859 


+53 


5.5(0.3) 


Fen 


13.63(0.05) 










Till 


11.34(0.10) 










Znn 


11.43(0.04) 










Sii 


14.88(0.12) 










Crn 


12.12(0.05) 










Pii 


13.11(0.04) 










Nin 


12.58(0.05) 










Ni 


< 15.16 










Mnii 


11.33(0.04) 










CI I 


12.90(0.03) 










Sin 


14.05(0.03) 










Mgii 


15.27(0.12) 










Sn 


14.05(0.10) 


3 


1.776538 


+18 


6.0(0.3) 


Fell 


13.46(0.04) 










Pn 


12.00(0.30) 










Znll 


11.81(0.02) 










Ni 


< 12.82 










Cm 


12.07(0.04) 










Mgn 


14.31(0.26) 










Niii 


12.27(0.10) 


6 


1.776964 


+64 


9.4(0.9) 


Fen 


13.23(0.07) 










Mnii 


11.53(0.02) 










Sin 


13.92(0.05) 










Sin 


14.40(0.01) 










Ni 


< 13.33 










Sn 


14.27(0.07) 






















Pn 


12.19(0.25) 






















Ni 


< 13.60 






















CI I 


12.50(0.06) 






















Mgn 


14.52(0.19) 














1 


1.776299 


-8 


15.0(1.0) 


Mgi 


11.72(0.07) 


4 


1.776684 


+34 


8.6(0.4) 


Mgi 


11.83(0.07) 


2 


1.776365 


-0.5 


6.0(0.2) 


Mgi 


12.08(0.03) 


5 


1.776814 


+48 


9.1(1.1) 


Mgi 


11.85(0.07) 










Ci 


13.12(0.02) 


6 


1.776947 


+62 


3.0(1.3) 


Mgi 


11.39(0.05) 


3 


1.776523 


-fl7 


3.4(0.3) 


Mgi 


11.77(0.01) 






















Ci 


12.72(0.02) 














Intermediate-ion transitions 


1 


1.776036 


-36 


8.4(0.4) 


Alni 


12.03(0.02) 


4 


1.776694 


+35 


15.1(1.4) 


Aim 


12.49(0.04) 


2 


1.776300 


-7 


8.1(1.1) 


Aim 


12.29(0.13) 


5 


1.776866 


+54 


7.2(0.5) 


Aim 


12.16(0.06) 


3 


1.776414 


+5 


7.8(1.4) 


Alni 


12.29(0.14) 















* Velocity relative to 2 = 1.776370 



ably this contamination. Finally, we obtain a more con- 
straining upper limit for A^(Ti+) from the Ti 11 A3384 line 
than the higher value deduced bv lProchaska et al] l)200l|) 
from the Ti 11 A1910 lines having a 3 times lower oscillator 
strength. 

The fitting solutions of the detected components in 
the low-ion metal-lines are presented in Table 0| A sur- 
prising particularity of this system known as a system 
exhibiting the largest dust depletion level of any DLA 
is that the fits of the refractory element lines require 6 
components (2 components at w ~ km s~^, see Fig. 0)), 
while the non-refractory and mildly refractory element 
lines are well fitted with 5 components only (1 component 
at w = km s~^, see Fig.0J. The refractory elements show- 
ing weak absorption lines, like the Mgii A1239,1240 and 



Cm A2056,2062 lines, are borderline cases, in the sense 
that the component 1 is detected but only marginally. 
Concerning the component 6, it is only observed in the 
strongest metal-line transitions. 

The detection of the Cli A1347 line i s exceptional. 
Cl° h as pre viously been obs e rved only by iLedoux et alJ 
l|2002|) and IProchaska et all l)2003l) in the DLA toward 
Q0551-366 and in the DLA toward FJ081240.6+320808, 
respectively. In the DLA system studied in this paper, we 
detect only the components 2 and 3 of the Cli A1347 line, 
the other components 4, 5 and 6 arc blended with Lya 
forest absorptions. However, since the line profiles of S 11 
A1250, Sill A1808 and Ci A1656 show very little absorp- 
tion at w > 25 km s~^ (only 20% of the total column 
density for Sii and Sill), we consider the Cli column den- 
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Fig. 4. Same as Fig.|21for the DLA toward Q1331+17. The zero velocity is fixed at z = 1.776370. 



sity derived from the components 2 and 3 as a value. We possible contribution from the blended components 4, 5 
adopt an error on iV(Cll) of 0.10 dex to account for the and 6. Nevertheless, the Cl° column density measurement 
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Fig. 3. Same as Fig. [T] but for Q1331+17. The verti- 
cal bars correspond to the wavelength centroids of the 
components used for the best fit, z = 1.776724 and 
z = 1.786345, referring to the DLA system and an ad- 
ditional absorber, respectively. The measured Hi column 
densities are logiV(Hi) = 21.14 ± 0.08 and 19.80 ± 0.10, 
respectively. 



cu 



o 

^ 3650 3700 3750 3800 

Observed Wavelength [A] 

Fig. 5. Same as Fig.[l]but for Q2231-00. The vertical bar 
corresponds to the wavelength centroid of the component 
used for the best fit left as a free parameter, z = 2.065774. 
The measured Hi column density is logiV(Hi) — 20.53 ± 
0.08. 



provides only a strict lower limit on the CI absolute abun- 
dance of [Cl/H] > —1.37, since the dominant state of CI 
should be Clii in DLAs (the ionization potential of Cl'^ 
is 13.01 eV which is lower than that of hydrogen). The 
CI II column density cannot unfortunately be determined 
in this DLA system due to the low rest-wavelengths of the 
CI II hues. 

The Ni column density is not easy to derive in this 
DLA. First, given the high H i column density of this DLA 
system, the Ni triplet at Aj-est ^ 1200 A is blended with 
the blue damping wing of the DLA Lya absorption line. 
We thus have to renormalize the spectral region cover- 
ing the Ni A1200 triplet with the fit of the Lya damp- 
ing wing profile, before fitting the Ni lines. Second, it 
is difficult to determine whether the derived N i column 
density is a lower or an upper limit. Indeed, the compo- 
nent 2 of all the detected Ni lines is saturated, even in 
the weakest Ni transition at Aiost = 1134.98 A. Thus, we 
are inclined to assume that the derived column density 
is a lower limit. But, when looking carefully at the Ni 
line profiles, they all show a stronger component 4 than 
the component 5, on the opposite to what is observed in 
other low-ionization transition profiles. This suggests that 
the Ni lines are likely blended with Lya forest absorp- 



tions. Therefore, we finally assume that the derived Ni 
column density is an upper limit. 

While fitting the Lya damping line, we found it nec- 
essary to include, besides the contribution of the DLA at 
z = 1.776724 corresponding to the redshift of the metal- 
line component 4, a second absorber shifted by about 1079 
km s^^ redwards of the DLA system (see Fig.|21l. The red- 
shift of this second absorber, z^hs = 1-786, is accurately 
determined thanks t o the detections of seve ral associated 
metal- lines (see also lLooez fc Ellisonll2003l) . The derived 
H I column densities of the DLA and the second absorber 
are log TV (Hi) = 21.14 ± 0.08 and 19.80 ± 0.10, respec- 
tively. They confirm the previous measurements made by 
IPettini et all l|l994j) . 



3.3. Q2231-00, Zabs = 2.066 

T his DLA systeiri was previously observed and analyzed 
bv lLu et al.rijl996|) a ndlfrochaska fc Wolfe ( 1999), and re- 
cently completed bv lProchaska et all |200lj) . Throughout 
their analysis they adopted the H i column density mea- 
sured by Lu & Wolfe (1994). We obtain here a new value 
for the logiV(Hi) of 20.53 ± 0.08 from the Lya damping 
line observed in the UVES spectra (see Fig. EJ. The fit 
was performed by fixing the 6- value at 20 km s^^ and by 
leaving the redshift as a free parameter. 

The UVES spectra provide column density measure- 
ments of the S+, P+, N°, Mn+ and Mg" ions (see Fig.E)), 
and confirm the iV(Si+), A^(Zn+), iV'(Cr+ ) and 7V(Ni+) 
measurements of lProchaska fc Wolfd ( 1999^ . The fit of the 
Mgi A2852 line clearly shows that the contamination of 
Znii A2026 by the Mgi A2026 profile is negligible. We note 
that we had to fit the Mgi profile with slightly different 
fitting parameters in comparison to the ones we have used 
for other low-ion metal lines (see Table ISJ. By combining 
the Fell A1608,1611 fines detected in the HIRES spec- 
tra and the Fell A2260,2344,2374 lines observed in the 
UVES spectra, we deduce a more accurate measurement 
of iV(Fe+). Finally, we obtain a revised value for iV(Tiii) 
of 12.66 ± 0.08 from the Tin A3067,3073 lines. We suc- 
ceed in deriving an upper limit on the column density of 
the intermediate-ion, Feiii, located in the Lya forest, by 
using the Aim lines observed in the HIRES spectra to 
constrain the fitting parameters (see Fig. and Table ISJ . 
The intermediate-ion lines show very similar profiles to 
the low-ion line profiles. 

The low-ion line profiles are characterized by a rela- 
tively complex velocity structure composed of 11 compo- 
nents presented in Table [3 and spread over 180 km s~^ 
in velocity space. But, only the 5 strongest components, 
the components 4, 5, 9, 10 and 11, are detected in the 
weak metal-lines. They contain only about 75% of the to- 
tal column density obtained by summing the contribu- 
tions of the 11 components. This shows that when com- 
puting the abundance ratios [X/Y], it is necessary to con- 
sider only the column densities of the velocity compo- 



10 M. Dessauges-Zavadsky et al.: A comprehensive set of elemental abundances in DLAs 



Table 5. Component structure of the Zabs — 2.066 DLA system toward Q2231— 00 



No 2abs 






Ion 


logAf(criogjv) 


No 


^ahs 






Ion 


logAf(o-iog]v) 






km s^^ 


km s^^ 










km s~^ 


km s~^ 






Low-ion transitions 




1 


2.064746 


-138 


6.6(0.5) 


Fell 


13.06(0.04) 


8 


2.065887 


-27 


4.9(0.1) 


Fen 


13.41(0.02) 










Sill 


13.53(0.12) 










Sin 


13.87(0.01) 


2 


2.064897 


-123 


9.0(0.7) 


Fell 


13.20(0.07) 










Sn 


13.80(0.26) 










Sill 


13.79(0.06) 










Ni 


13.59(0.22) 


3 


2.065087 


-105 


13.8(1.4) 


Fen 


13.58(0.04) 


9 


2.066072 


-9 


5.5(0.4) 


Fen 


14.00(0.05) 










Sin 


13.97(0.10) 










Znn 


11.71(0.06) 


4 


2.065266 


-88 


8.6(0.2) 


Fell 


14.19(0.01) 










Crn 


12.12(0.08) 










Znn 


11.63(0.04) 










Niii 


12.73(0.03) 










Cm 


12.46(0.03) 










Mnll 


11.92(0.05) 










Niii 


13.04(0.10) 










Sin 


14.59(0.04) 










Mnii 


11.97(0.03) 










Tin 


11.88(0.18) 










Sin 


14.57(0.04) 










Sn 


14.33(0.13) 










Sii 


14.19(0.13) 










Pn 


< 12.81 










Pn 


< 12.90 










Ni 


14.43(0.06) 










Ni 


13.95(0.11) 


10 


2.066161 





3.8(0.2) 


Fen 


14.26(0.03) 


5 


2.065422 


-72 


4.2(0.5) 


Fen 


13.59(0.03) 










Znn 


11.98(0.03) 










Znn 


10.92(0.17) 










Crn 


12.60(0.02) 










Crii 


11.77(0.13) 










Niii 


13.01(0.02) 










Niii 


12.44(0.17) 










Mnll 


12.16(0.03) 










Mnii 


11.46(0.07) 










Sin 


14.83(0.03) 










Sin 


13.90(0.06) 










Tin 


12.58(0.06) 










Sn 


13.87(0.16) 










Sn 


14.66(0.13) 










Pn 


< 12.32 










Pn 


< 13.03 










Ni 


13.22(0.03) 










Ni 


14.76(0.22) 


6 


2.065546 


-60 


6.4(0.2) 


Fen 


13.58(0.02) 


11 


2.066275 


-fll 


5.4(0.1) 


Fen 


13.72(0.01) 










Sin 


13.83(0.08) 










Crn 


12.07(0.08) 










Sn 


13.99(0.14) 










Niii 


12.70(0.04) 










Ni 


13.25(0.13) 










Mnll 


11.60(0.06) 


7 


2.065698 


-45 


4.7(0.3) 


Fen 


13.18(0.02) 










Sin 


13.94(0.03) 










Sin 


13.57(0.01) 










Sn 


14.06(0.21) 










Sn 


13.89(0.25) 










Pn 


< 12.74 










Ni 


13.30(0.23) 










Ni 


13.50(0.23) 


1 


2.064788 


-134 


7.6(3.0) 


Mgi 


11.04(0.12) 


6 


2.065903 


-25 


7.0(0.7) 


Mgi 


11.36(0.03) 


2 


2.064981 


-115 


7.9(2.1) 


Mgi 


11.12(0.09) 


7 


2.066064 


-9 


5.2(0.4) 


Mgi 


11.72(0.05) 


3 


2.065233 


-91 


8.8(0.5) 


Mgi 


11.87(0.03) 


8 


2.066155 


-0.6 


3.1(1.3) 


Mgi 


12.11(0.06) 


4 


2.065461 


-68 


14.7(1.7) 


Mgi 


11.65(0.05) 


9 


2.066274 


-fll 


4.6(1.0) 


Mgi 


11.59(0.06) 


5 


2.065692 


-46 


3.0(1.0) 


Mgi 


11.05(0.06) 














Intermediate-ion transitions 


1 


2.064758 


-137 


3.0(1.2) 


Aim 


11.72(0.07) 


6 


2.065648 


-50 


10.7(2.3) 


Aim 


11.88(0.07) 










Fein 


< 13.64 










Fein 


< 13.83 


2 


2.064909 


-122 


6.6(1.5) 


Aim 


12.00(0.05) 


7 


2.065925 


-23 


8.5(2.7) 


Aim 


12.10(0.10) 










Fein 


< 13.92 










Fein 


< 13.83 


3 


2.065123 


-101 


8.5(2.3) 


Aim 


12.01(0.10) 


8 


2.066076 


-8 


3.3(1.0) 


Aim 


12.04(0.19) 










Fein 


< 13.78 










Fein 


< 13.59 


4 


2.065274 


-87 


6.5(0.7) 


Aim 


12.50(0.03) 


9 


2.066165 





5.1(2.5) 


Aim 


12.51(0.09) 










Fein 


< 13.59 










Fein 


< 13.74 


5 


2.065435 


-71 


3.0(1.6) 


Aim 


11.97(0.02) 


10 


2.066290 


+13 


4.8(1.3) 


Aim 


11.96(0.13) 










Fein 


< 13.77 










Fein 


< 13.59 



* Velocity relative to z = 2.066161 



nents detected in both the X and Y profiles to avoid un- 
der /overestimations of the relative abundances. 

An additional difficulty which appears in absorption 
systems with complex metal-line profiles is that the prob- 
abihty of blending in the Lya forest is higher over a larger 



velocity interval. For this reason it is more sensible to con- 
sider the measured P+ and N° column densities as upper 
limits, the Pil and Ni lines being located in a region of 
the Lya forest with numerous absorption lines (see Fig.EJ. 
The S"^ column density is a borderline case, and unfortu- 
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Same as Fig.Hlfor the DLA toward Q2231-00. The zero velocity is fixed at z ^ 2.066161. 
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nately only the Sii A1253 line from the S+ triplet can be ily blended with Lya clouds. We consider the measured 
used to determine iV(S+), the two other lines are heav- 
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Table 6. Component structure of the Zabs — 2.431 DLA system toward Q2343+12 



No 






biat) 


Ion 


logiV(aiogjv 






km s^^ 


km s^^ 






Low-ion transitions 


1 


2.428181 


-154 


3.0(1.0) 


Fell 


12.44(0.03) 










Oi 


14.38(0.03) 










Sin 


13.59(0.02) 


2 


2.428277 


-146 


3.0(0.9) 


Fell 


12.94(0.02) 










Oi 


14.98(0.04) 










Sin 


13.81(0.06) 


3 


2.428383 


-136 


4.1(0.3) 


Fell 


13.18(0.01) 










Oi 


14.52(0.02) 










Sin 


14.12(0.06) 


4 


2.428499 


-126 


4.4(0.5) 


Fell 


12.49(0.04) 










Oi 


13.30(0.19) 










Sin 


13.37(0.01) 


5 


2.428749 


-104 


8.5(0.3) 


Fell 


12.74(0.01) 










Oi 


14.51(0.02) 










Sin 


13.08(0.05) 


6 


2.428884 


-93 


3.0(1.0) 


Fell 


12.68(0.02) 










Oi 


14.29(0.03) 










Sin 


12.95(0.02) 


7 


2.429009 


-82 


5.9(0.4) 


Fell 


12.45(0.02) 










Oi 


14.16(0.03) 










Sin 


12.87(0.02) 


8 


2.429249 


-61 


5.1(0.4) 


Fell 


12.50(0.02) 










Oi 


14.15(0.03) 










Sin 


13.14(0.05) 


9 


2.429355 


-51 


3.0(0.8) 


Fell 


12.85(0.03) 










Oi 


14.54(0.02) 










Sin 


13.71(0.02) 


10 


2.429500 


-39 


5.2(0.2) 


Fell 


12.65(0.01) 










Oi 


14.22(0.03) 










Sin 


13.28(0.09) 


11 


2.429702 


-21 


4.4(0.5) 


Fell 


12.20(0.02) 










Oi 


13.97(0.05) 










Sin 


12.81(0.06) 


12 


2.429942 





7.3(0.2) 


Fen 


13.28(0.01) 










Oi 


14.89(0.01) 










Sin 


13.73(0.11) 


13 


2.430086 


+13 


9.6(0.8) 


Fen 


12.86(0.06) 










Oi 


14.59(0.02) 










Sin 


13.33(0.01) 


14 


2.430394 


+40 


9.7(0.2) 


Fen 


13.19(0.01) 










Oi 


14.75(0.01) 










Sin 


13.68(0.01) 



* Velocity relative to z ^ 2.429942 



No 2abs -Urei b{at) Ion logN{aiogN) 

km s~^ km s~^ 



15 


2.430674 


+64 


2.8(0.2) 


Fen 


12.38(0.03) 










Oi 


13.93(0.05) 










Sill 


13.13(0.01) 


16 


2.430808 


+76 


3.2(0.1) 


Fen 


12.65(0.02) 










Oi 


14.37(0.03) 










Sin 


13.32(0.01) 


17 


2.431018 


+94 


8.3(0.4) 


Fen 


13.45(0.02) 










Oi 


15.27(0.05) 










Sin 


14.01(0.09) 


18 


2.431157 


+106 


5.1(0.4) 


Fen 


13.98(0.04) 










Znii 


11.72(0.03) 










Crii 


12.37(0.03) 










Mnll 


11.82(0.04) 










Oi 


15.76(0.06) 










Sin 


14.55(0.03) 










An 


12.85(0.18) 










Ni 


14.31(0.22) 










Pii 


< 12.67 


19 


2.431288 


+118 


6.2(0.6) 


Fen 


14.09(0.03) 










Znii 


11.87(0.02) 










Crii 


12.38(0.03) 










Mnll 


11.89(0.03) 










Oi 


15.90(0.05) 










Sin 


14.71(0.02) 










An 


12.86(0.10) 










Ni 


14.32(0.15) 










Pii 


< 12.55 


20 


2.431439 


+131 


4.7(0.4) 


Fen 


13.82(0.03) 










Znii 


11.67(0.03) 










Crii 


12.20(0.05) 










Mnll 


11.64(0.05) 










Oi 


15.80(0.05) 










Sin 


14.54(0.03) 










An 


12.44(0.22) 










Ni 


13.58(0.25) 










Pii 


< 12.45 


21 


2.431570 


+142 


5.7(0.2) 


Fen 


13.48(0.04) 










Oi 


15.09(0.07) 










Sin 


13.76(0.16) 


22 


2.431826 


+165 


6.9(0.1) 


Fen 


13.02(0.01) 










Sin 


13.77(0.01) 


23 


2.432126 


+191 


6.9(0.5) 


Fen 


12.27(0.03) 










Sin 


12.58(0.03) 



7V(S+) as a value, assuming that the adopted error on 
iV(S+) covers the possible Hi blendings. 



3.4. Q2343+12, Zabs = 2.431 



This D LA system has first been studied bv lSargent et ahl 
l)l988|) . The chemical abundance measurements obtained 
from HIRES spectra were r nainly used in s tatistical 
samples (e.g. ILu et all Il998l IProchaska ATWolfe. .1998c 



Vladilo et a,1.l2nfii Ti.ecentlv.ln'Odorico et all and 



Vladilo et al.1 l(2003l) reported some additional abundance 



measurements from UVES spectra. We present here the 
first complete set of elemental abundances of this DLA. 

From our UVES spectra, we measured the column den- 
sities of Oi, An, Ni, Feii, Znii, Cm and Mnii, and ob- 
tained an upper limit on the column density of P ii con- 
taminated by the Lya clouds. The metal absorption pro- 
files are characterized by an extremely complex velocity 
structure extended over 350 km s~^ in velocity space and 
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Table 7. Component structure of the Zabs — 2.431 DLA system toward Q2343+12 — Continued 



No 


^ahs 


km s^^ 


km s~^ 


Ion 


logiV((Jiogjv) 


No 


^ahs 


"rcl 

km s~^ 


km s~^ 


Ion 


log7V(criog]v) 


Intermediate-ion transitions 


1 


2.428178 


-154 


6.5(2.2) 


Aim 
Fem 
Nil 
Siii 


11.93(0.25) 

13.06(0.06) 
13.99(0.01) 
< 14.17 


7 


2.430392 


+39 


10.8(0.5) 


Aim 
Fem 
Nil 
Siii 


11.92(0.03) 

< 13.40 
13.66(0.02) 

< 13.84 


2 


2.428288 


-144 


3.6(1.6) 


Aim 
Fem 
Nil 


11.90(0.27) 

13.19(0.04) 
14.48(0.21) 


8 


2.430696 


+66 


10.8(0.4) 


Aim 

Nil 

Siii 


11.72(0.05) 

13.71(0.01) 
< 13.90 










Siii 


< 13.56 


9 


2.431065 


+98 


15.0(0.2) 


Aim 


12.08(0.26) 


3 


2.428426 


-133 


9.3(0.6) 


Aim 
Fem 


12.59(0.03) 
13.54(0.02) 










Nil 
Siii 


14.04(0.01) 
< 14.35 










Nil 


14.36(0.01) 


10 


2.431185 


+109 


8.9(1.0) 


Aim 


12.45(0.08) 










Siii 


< 14.22 










Nil 


> 14.25 


4 


2.429358 


-51 


8.3(0.5) 


Aim 


11.51(0.07) 










Siii 


< 13.71 










Fem 


< 10.00 


11 


2.431323 


+121 


6.7(0.9) 


Aim 


12.39(0.08) 










Nil 


13.89(0.04) 










Nil 


> 14.11 










Siii 


< 13.97 










Siii 


< 14.36 


5 


2.429487 


-40 


19.5(1.3) 


Aim 
Fem 
Nil 


11.83(0.05) 

< 14.09 
13.80(0.04) 


12 


2.431443 


+131 


3.0(1.5) 


Aim 

Nil 

Siii 


11.92(0.08) 

> 13.66 
< 13.45 










Siii 


< 13.62 


13 


2.431545 


+140 


10.5(1.5) 


Aim 


12.23(0.06) 


6 


2.429995 


+5 


14.4(0.8) 


Aim 
Fem 


12.09(0.02) 
< 13.36 










Nil 
Siii 


> 14.27 
< 14.59 










Nil 


13.71(0.02) 


14 


2.431849 


+167 


7.7(0.2) 


Aim 


11.72(0.04) 










Siii 


< 13.79 


15 


2.432135 


+192 


15.0(2.1) 


Nil 
Siii 
Nil 


14.01(0.01) 

< 14.34 
13.72(0.01) 



* Velocity relative to z = 2.429942 



X in 



3 



g 




^ 3505 3510 3515 3520 3525 3530 
Observed Wavelength [A] 

Fig. 7. Same as Fig. [T] but for Q2343+12. The Voigt 
fit was performed w i th the Hi column density measured 
by iD'Odorico et all l|20n2l Sect. 3.2) from the Lya line 
(logiV(Hi) = 20.35 ± 0.05). The vertical bar corresponds 
to the wavelength centroid of the component used for the 
best fit by these authors, z — 2.43125. 



composed of 23 components (see Fig.|Sland TableE)). The 
dominant feature is the group of 5 components around 
V ~ +120 km s^^, with the components 18, 19 and 20 
being the strongest ones. These components are heavily 
saturated in strong metal-lines, but are the only compo- 
nents detected in weak metal-lines. They contain about 
62% of the total column density. 



The Lya line is outside our UVES wavelength 
coverage, therefore we adopt the H i column density 
measurement (log NjH i) = 20.35 ± 0.05) obtained by 
ID'Odorico et all l)2002l) . Fig. Q shows the fit of the Ly/3 
line which confirms this iV(H*^) measurement. This rela- 
tively low H'' column density and the presence of a strong 
N II A1083 line as well as of other intermediate-ion lines — 
Aim A1854, Fem A1122 and Siii A1012 - with very sim- 
ilar profiles to the low-ion metal line profiles (see Fig. [S] 
and Table [Tj) suggest that some photoionization correc- 
tions may be required in this DLA system. A discussion 
on the photoionization corrections is presented in Sect. 01 

4. Photoionization corrections 

The photoionization effects have to be carefully examined 
in absorption metal-line systems, when the main goal is 
to provide a complete and unbiased interpretation of the 
chemical abundance patterns. Indeed, as we are studying 
gas-phase abundances, a fraction of the gas may be ion- 
ized and we need to determine this fraction to be sure 
that we are measuring the intrinsic abundances. The dust 
depletion effects also affect the gas-phase abundance mea- 
surements. The manner we tackle this issue is described 
in the next Section. 
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Fig. 8. Same as Fig.Elfor the DLA toward Q2343+12. The zero velocity is fixed at z = 2.429942. 
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It is generally assumed that in the DLA systems the tion state in H i regions is the neutral one for elements 
ionization fraction is low, so that the dominant ioniza- with the first ionization potential > 13.6 eV (e.g. O", N") 
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and the singly ionized one for elements with the first ion- 
ization potential < 13.6 eV and the second > 13.6 eV 
(e.g. Fe+, Si+). The reason for this is that the bulk of 
the H° gas is self-shielded from hv > 13.6 eV photons, 
but transparent to hv < 13.6 eV photons. Under this as- 
sumption, the correction factors for ionization are small, 
and thus the column densities of low-ions are equal to 
the column densities of elements, e.g. Af(Si+) ~ iV(Si). 
However, the detection of intermediate-ionization transi- 
tions, like Al^"*", in the DLA systems provided doubts on 
the ionization levels in DLAs. The ionization potential of 
A1+ is 18.8 eV, i.e. greater than that of hydrogen, therefore 
A1++ is likely present in ionized and not in neutral gas, 
since photons with hv > 18.8 eV cannot so easily pen- 
etrate gas clouds with large H i column densities. In ad- 
dition, the observations show a good correlation between 
the velocity stru ctures of Al ill and singly ionized and neu- 
tral species (e.g.l Lu et al.' ll996t IProchaska fc Wolfe' 1999; 
IWolfe fc Prochaska .2000) . To explain the similarity be- 
tween the Aim and low -ion ization species line pro files, 
iHowk fc SembachI (|l299) and'izot ov fc ThuanI ^999) pro- 
posed that these lines originate in the same ionized region 
or in a mix of neutral and ionized clouds, and stressed the 
importance of abundance corrections for ionization effects. 
Furthermore, indirect arguments to consider ionized gas 
in DLAs came from studies of the warm ioni zed medium 
in th e Milk y Way and other nearby gal axies (H owk et al.l 
Il999e .Sembach et allEoOol: I Jenkins et al 20001 . 

Several authors have investigated the ionization effects 
in the DLA systeins llv'iegasll99.4lHowk fc Sembachll999t 
llzotov et allboOli IVladilo et al.ll200ll) . Although the dif- 
ferent approaches used to deal with the problem of pho- 
toionization in DLAs have led to slightly diff'erent conclu- 
sions, it was generally found that ionization corrections in 
DLAs are negligible, being below the measurement errors, 
except for some particul ar systems like the DLA toward 
GB 1759-^7539 studied bv lProchaska et alJ ll2002bl) . While 
the authors stress that this particular DLA system has 
several characteristics which separate it from the major- 
ity of DLAs, its properties highlight the importance of as- 
sessing the ionization state of each DLA system. The pho- 
toionization model computations are generally performed 
using the CLOUDY photoionization equilibrium software 
package (e.g. iFter land ct al. 1998). To avoid computin 
heavy photoionization models, IProchaska et alJ l|2002 
defined a number of photoionization diagnostics which 
provides a qualitative "first-look" analysis of the ioniza- 
tion state in a DLA system and of the level of required 
ionization corrections. In the following sub-Sections, we 
apply these photoionization diagnostics to each of the four 
DLAs studied and qualitatively evaluate the importance 
of ionization corrections. 



4.1. QOlOO+13, Zabs = 2.309 

According to the anti-correlation between logA^(H'') 
and logAr(Al++)/7V(Al+) found bv iVla,dilo et al... (.200 1,V 



where the A1++/A1+ column density ratio is assumed to 
be at a first approximation an indicator of the ionization 
level in the gas, the high Hi column density, logiV(Hi) 

— 21.37 ± 0.08, of this DLA system suggests that the ion- 
ization corrections are low. The following three additional 
and more reliable ionization indicators are at our disposal 
in this DLA system. 

Thanks to the high efficiency of UVES in the blue, 
we have the spectral coverage for this DLA system of the 
rarely detected intermediate-ion transitions, Fe ill and N ii 
at Arost — 1122 and 1083 A, respectively. The column 
density ratios Fc++/Fc+ and are the most reli- 

able ionization indicators, since the recombination coef- 
ficient of Al"*" is rather imcertain, being likely overesti- 
mated ijNussbaumer fc Storevlll98(^ . The observed Feiii 
and N ii lines show the same profiles as the low-ion profiles 
and are well fitted with the fitting parameters deduced 
for the low-ion lines. Located in the Lya forest, one of 
the 2 components is clearly blended in both intermediate- 
ion lines (see Fig. [21), therefore we derived only upper 
limits on their column densities. These limits lead to 
logiV(Fe++)/iV(Fe+) < -1.54 and logN{N+)/ N(N°) < 

— 1.24. Acco rding to the ph otoionization diagnostics es- 
tablished by IProchaska et al.l l|2002blh . in a system where 
Fe++/Fe+ < -1.60 dex and N+/N0 < -1.00 dex, the 
ionization fraction, x, defined as the ratio of H''" over 
(H°-|-H+), is expected to be lower than 10%. Thus, the 
measured Fe and N"*" column densities are in agree- 
ment with the expectations for a DLA with a H i column 
density higher than 10^^ cm~^. 

Another ionization indicator that we can consider in 
this system is argon. Indeed, Ar has a very high photoion- 
ization cross-section for photons with energy higher than 
13.6 eV which is about ten times higher than the one 
of Hi, and thus Ar is very sensitive to ionization effects 
( Sofia fc ,Ienkinslll998f ). The comparison of the Ar abun- 
dance with the abundance of another a-element — O, S, 
or Si — provides indication on the ionization level in the 
gas. In this s ystem, we measur e [Ar/S ] = — 0.20±0.14, and 
according to IProchaska et alJ l)2002b|) . an observed value 
of [Ar7S+] > -0.20 dex indicates x < 10%. 

In summary, the high H i column density, the 
A1++/A1+, Fe++/Fe+ and N+/N° column density ratios, 
and the Ar'^/S"'" abundance ratio, all point to a low ioniza- 
tion fraction in this DLA system, and thus to low ioniza- 
tion corrections. The ionization corrections should indeed 
be lower than 0.1 dex for all the elements, except probably 
for Ar. 

4.2. Q1331+17, zabs = 1-776 

This is another very high H i column density system 
with log7V(Hi) = 21.14 ± 0.08. The high Hi column 
densit y suggests that th e ionization level is low in this 
DLA llVladilo et alJl200ll) . However, we do not detect any 
very reliable ionization indicator in this system to con- 
firm this statement. Only the A1++/A1+ column density 
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ratio is available for which we derive an upper limit of 
< —0.76 dex from the HIRES spectra. This limit indi- 
cates that the ionization corrections should be low in this 
system. 

4.3. Q2231-00, z^bs = 2.066 

Two ionization indicators are accessible in this DLA with 
a Hi column density of logA^(Hi) = 20.53 ± 0.08. The 
first one is the A1++/A1+ column density ratio obtained 
from the HIRES spectra, and the second one is the up- 
per limit on the Fe++/Fc+ column density ratio obtained 
from the UVES spectra. However, the Alii line is so heav- 
ily saturated that no reliable A^(A1"'') lower limit can be 
deduced. We thus derive an estimation of the A1++/A1+ 
column density ratio by using Si"*" as a proxy of A1+ assum- 
in g the logA^(Si ~ '") ver sus log7V(Al+) correlation found 
hv lviadilo et"an l)200lf) . We obtain logiV(Al++)/iV(Al+) 
= —0.76 ± 0.19. With regard to the second ionization in- 
dicator, the Feiii A1122 line is located in the Lya forest 
and is strongly blended with Hi clouds. Consequently, the 
deduced Fe"''^ column density upper limit is not stringent. 
We get \ogN(Fe++)/N(Fe+) <C -0.09. 

Unfortunately the Ar i lines as well as the N ii line are 
beyond the quasar flux cut-off, no other ionization indica- 
tor is hence accessible in this DLA system. If we trust the 
derived A1++/A1+ column density ratio measurement, the 
ionization corrections should not be significant. 

4.4. Q2343+12, z^bs = 2.431 

The low Hi column density of logiV(Hi) = 20.35 ± 0.05 
suggests that the ionization corrections might be relevant 
in this system. Several strong intermediate-ion lines with 
very similar profiles to the low-ion line profiles are ob- 
served in this system, namely Aim, Fciii, Nil and Siii 
(see Fig. |Hll. This is indeed a first indication that pho- 
toionization is playing an important role here. The second 
straightforward indication is provided by the Ar/0 and 
Ar/Si abundance ratio measurements. They are very low 
in this system. We find [Ar/0] = -0.72 ±0.17 and [Ar/Si] 
= —0.81 lb 0. 1 6. in agreement with the measurements by 
IVIadilo et alJ ll200,'^ . The Ar/S i abundance ratio in this 
DLA system is even lower than the [Ar/Si] abundance ra- 
tio of —0.68 ± 0.0 4 measured in the DLA system toward 
GB 1759-^7539 bv lProchaska et all l)2002b|) whose analysis 
showed important ionization corrections. 

The measured logiV(Al++)/7V(Al+) = -0.62 ± 0.19 
obtained by taking Si+ as a pr oxy of A1+ accord- 
ing to the correlation identified by 'Vladilo et alj lj200ll) , 
logiV(Fe++)/iV(Fe+) < 0.37 derived from the Feili A1122 
line which is only partially blended with Lya clouds (in the 
components 4 and 5), and log A^(N+)/A^(N°) > -O.lOcon- 
firm that the ionization corrections are important in this 
DLA system. Indeed, a system with Fe++/Fe+ > —1 dex 
and N+/N° > —0.2 dex is expect ed to have an ionization 
fraction, x, higher than 50% IjProchaska et al...2002b.) . We 



suspect that the ionization corrections in this DLA system 
are of similar magnitude as the ones of the DLA toward 
GB 1759+7539, namely between 0.1 to > 0.5 dex depend- 
ing on the ion. Such high ionization corrections prevent us 
from deriving reliable chemical abundances in this DLA 
system. Therefore, we do not include this DLA in the 
following Sections where we discuss the chemical abun- 
dances of individual systems and make a detailed com- 
parison with chemical evolution models. 

5. Chemical abundances 

The small sample of DLAs studied in this paper is unique, 
because by combining the UVES-VLT data with the 
HIRES-Keck data we could measure the abundances of 
up to 15 elements - N, O, Mg, Al, Si, P, S, CI, Ar, Ti, 
Cr, Mn, Fe, Ni, Zn. This contrasts with the majority of 
DLAs for which only a handful of element s (Si, Fe. occa- 
sionally Cr, Zn, Ni) is usually detected fe.g.lLu et al.ll996l: 
IProchaska, WoifgiT99ilp7^cha,ska et a,lJ200lh. This low 
amount of information on individual systems has until now 
severely limited the interpretation of the DLA abundance 
patterns. 

5.1. Dust content 

The interpretation of the elemental abundance patterns in 
DLAs is by far not straightforward, the principal difficulty 
is to disentangle the nucleosynthetic contributions from 
the dust depletion effects. Several pieces o f evidence show 
that dust is indeed present in DLAs (e.s. IPei et al.lll99lt 
IPettini et a1.l [T99i Ihou et alJ EoOll IProcha,ska Wolfel 
'20021. Therefore, as we are measuring gas-phase elemen- 
tal abundances in DLAs, in presence of dust the observed 
abundances may not represent the intrinsic chemical com- 
position of the system if part of the elements is removed 
from the gas to the solid phase, as it happens in the inter - 
stellar medium of our Galaxy ijSavage fc Sembachlll996|) . 
Consequently, the refractory elements (e.g. Si, Fe, Cr, Ni) 
preferentially incorporated into dust grains are not the 
best diagnostic elements and their relative ratios have to 
be cautiously interpreted, since their differential depletion 
can mimic the expected nucleosynthetic abundance pat- 
terns. 

In Sect. 2] we have seen that at least the ionization 
effects are playing a negligible role in three out of the four 
DLAs studied in this paper. To constrain now the dust 
depletion effects, the access to a large number of elements 
is crucial, since elements with the same nucleosynthetic 
origin provide information on the dust depletion level in 
DLAs. 

Fe-peak elements and Zn. When comparing the absolute 
abundances of different iron-peak elements, Cr, Fe, Ni, and 
of Zn measured in the DLAs toward QOlOO+13, Q1331-bl7 
and Q2231— 00, wc notice variations, while in the Galactic 
stars, Cr, Fe, Ni and Zn track each other and have solar 
values relative to Fe (apart from Zn which seems to be 
enhanced by at maximum 0.1 — 0.2 dex relative to Fe, 
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Table 8. Summary of the absolute abundances in the three DLA systems studied 



Quasar 


QOlOO+13 


Q1331+17 


Q2231-00 




2.309 


1.776 


2.066 


logAf(Hi) 


21.37(0.08) 


21.14(0.08) 


20.53(0.08) 



logAf(Znii) 


Zn/H]obs 


12.47(0.01) 


-1.57(0.09) 


12.54(0.02) 


-1.27(0.09) 


12.30(0.05) 


-0.90(0.10) 


log7V(Fell) 


Fe/H 
Fe/H' 
Fe/H 
Fe/H' 

fc/h; 


obs 

cor (EDO) 
cor Ell) 

cor (SOO) 

cor (Sll) 


15.09(0.01) 


-1.78(0.08) 
-1.63(0.10) 
-1.63(0.10) 
-1.52(0.10) 
-1.52(0.10) 


14.63(0.03) 


-2.01(0.09) 
-1.23(0.12) 
-1.15(0.13) 
-0.95(0.13) 
—0.95(0.13) 


14.83(0.03) 


-1.20(0.09) 
-0.83(0.14) 
-0.83(0.14) 
-0.69(0.14) 
-0.69(0.14) 


log7V(Crll) 


Cr/H 
Cr/H 
Cr/H 
Cr/H 
Cr/H 


obs 

cor (EOO) 
cor (Ell) 
cor SOO) 
cor (Sll) 


13.37(0.01) 


-1.69(0.08) 
-1.61(0.09) 
-1.60(0.10) 
-1.53(0.10) 
— 1 51 (0 10) 


12.95(0.03) 


-1.88(0.09) 
-1.23(0.13) 
-1.14(0.15) 
-1.07(0.12) 
—0 94(0 1 3) 


13.00(0.04) 


-1.22(0.09) 
-0.94(0.13) 
-0.94(0.13) 
-0.80(0.13) 
—0 82(0 1 3) 


logA''(Nill) 


Ni/H 
Ni/H 
Ni/H 
Ni/H 
Ni/H 


obs 

cor (EOO) 
cor (Ell) 
cor (SOO) 
cor (Sll) 


13.87(0.01)" 


-1.75(0.08)" 

-1.62(0.10) 

-1.60(0.10) 

-1.52(0.10) 

-1.50(0.10) 


13.44(0.08) 


-1.95(0.11) 
-1.23(0.15) 
-1.09(0.17) 
-0.95(0.14) 
-0.88(0.15) 


13.54(0.06) 


— 1.24(0.10) 
-0.91(0.14) 
-0.87(0.15) 
-0.78(0.14) 
-0.73(0.14) 


loo- N(Mn u) 


Mn/H]obet 






12 50(0 03) 


— 2 1 7(0 09) 


1 2 59(0 04) 


— 1 47(0 OQ) 


loB TVfSiii) 


Si/H] 
Si/H 
Si/H 
Si/H] 
Si/H] 


3bs 

.or (EOO) 
3or (Ell) 
=or (SOO) 
.or (Sll) 


> 14.72" 


> —2.21" 

> -2.21 

> -2.21 

> -2.21 

> —2.21 


15 30(0 01) 


— 1 40(0 08) 
-1.22(0.11) 
-1.16(0.12) 
-0.98(0.11) 

— 1 01 (0 1 1 "1 


15 29(0 04) 


—0 80(0 W) 
-0.78(0.10) 
-0.78(0.10) 
-0.74(0.11) 
—0 74(0 1 Ol 


Inp- /V(S it) 
log i\ ii; 


S/H]obs 


1 5 nqfo 06) 


— 1 48(0 1 1 ) 


1 5 ns(o 1 1 ) 


—1 26(0 ^4) 


1 5 1 0(0 1 5) 


—0 63(0 1 7) 


log7V(Mgii) 


Mg/B 
Mg/H 
Mg/H 
Mg/H 
Mg/H 


obs 
dcor (EOO) 
[]cor (Ell) 
dcor SOO) 
dcor (Sll) 


15.57(0.09) 


-1.38(0.11) 
-1.38(0.11) 
-1.38(0.11) 
-1.37(0.11) 
-1.37(0.11) 


15.53(0.14) 


-1.19(0.15) 
-1.05(0.16) 
-0.90(0.19) 
-0.87(0.17) 
-0.76(0.18) 






logAr(Tiii) 


Ti/H] 


* 

obs 


< 12. 21** 


< -2.10*" 


< 11.34 


< -2.74 


12.66(0.08) 


-0.81(0.11) 


logAf(Ari) 


Ar/H]ob= 


14.21(0.12) 


-1.68(0.15) 










logAr(Ni) 


N/H]obs 


15.03(0.10) 


-2.31(0.14) 


< 15.23 


< -1.88 


< 15.02 


< -1.48 


logiV(Pii) 


P/H]ob, 


13.05(0.09) 


-1.85(0.13) 


13.25(0.10) 


-1.42(0.13) 


< 13.51 


< -0.55 


logAf(Alii) 


Al/H] 


obs 






> 13.74" 


> -1.89" 






logiV(Cli) 


Cl/H] 


obs 






13.05(0.10) 


> -1.37^= 







"^rochaska & Wolfe (1999) 
'' IProch aska ct alj (i2001i) 



The [Cl/H] absolute abundance derived from the CI I column density has to be considered as a strict lower limit, since the 
dominant state of CI should be Clll in DLAs. 

^ We do not compute the dust corrections for [Mn/H], since this element when analyzed is compared with Fe, and these two 
elements have very similar dust depletions. 

* We do not compute the dust corrections for [Ti/H], since this element when analyzed is compared with Fe, and in this 
case dust depletion and nucleosynthesis tend to work in the opposite sense (see Paper I). 

In the dust correction models EOO and Ell, Vladilo j 2002a) assumes that the intrinsic [Zn/Fe] ratio is equal to — +0.10 dex. 
In the dust correction models SOO and Sll. lVladild ll2002aD assumes that the intrinsic [Zn/Fe] ratio is equal to = +0.00 dex. 



see IProchaska et aLllioont iMishenina e t al."2002'). These 
absolute abundance variations are in line with the vari- 
ations observed in t he Galactic ISM due to the differ- 
ential depletion (e.g. ISavaee fc Semba' ch''l996V and are 
thus suggestive as being the result of a differential deple- 
tion present in the DLA systems. Hence, the measured 
[Zn/Fe, Cr,Ni] ratios, when assumed as being the result of 
the different degree of incorporation into dust grains of Zn, 
a volatile element (almost undepleted), and of Fe, Cr and 
Ni, refractory elements, clearly show that dust depletion 



has to be seriously taken into account in the DLA sys- 
tems toward Q1331+17 and Q2231— 00, for which we find 
[Zn/Fe] = +0.75±0.05 and +0.45±0.07, respectively. The 
DLA toward Q1331+17, in particular, exhibits one of the 
largest dust depletion level of any DLA. On the contrary, 
the DLA toward QOlOO+13 shows a low [Zn/Fe] ratio of 
+0.25 ±0.04, and thus requires negligible dust corrections. 

Alpha- elements. Another similar way to highlight the pres- 
ence of dust in the DLA systems is the comparison of 
the absolute abundances of different a-elements. In the 
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Table 9. Summary of the relative abundances in the three DLA systems studied 



Quasar 


QOlOO+13 


Q1331+17 


(j2231— 00 




2.309 


1.776 


2.066 


logAf(Hi) 


21.37(0.08) 


21.14(0.08) 


20.53(0.08) 




Zn/Fe 


obs 


+0.25(0.04) 


+0.75(0.05) 


+0.45(0.07) 




Zn/Fe 


]cor (E) 


+0.10 


+0.10 


+0.10 




Zn/Fe 


]cor (S) 


+0.00 


+0.00 


+0.00 




Ni/Fe 


obs 


+0.07(0.02) 


+0.07(0.09) 


+0.07(0.07) 




Ni/Fe 


cor (EOO) 


+0.05(0.02) 


+0.01(0.02) 


+0.03(0.04) 




Ni/Fe 


CO. (Ell) 


+0.07(0.02) 


+0.07(0.09) 


+0.07(0.07) 




Ni/Fe 


CO. (SOO) 


+0.04(0.01) 


+0.00(0.01) 


+0.02(0.03) 




'Ni/Fe 


cor (Sll) 


+0.07(0.02) 


+0.08(0.09) 


+0.07(0.07) 




Mn/F 


e]obs 




-0.15(0.04) 


-0.16(0.05) 




Si/Fe 


obs 


> -0.43 


+0.61(0.03) 


+0.40(0.05) 




Si/Fe 


cor (EOO) 


> -0.58 


+0.01(0.03) 


+0.05(0.09) 




Si/Fe 


cor (i^ll) 


> -0.58 


-0.01(0.05) 


+0.05(0.09) 




Si/Fe 


cor (bUUj 


> -0.67 


-0.02(0.01) 


-0.05(0.07) 




Si/Fe 


cor (bll) 


> -0.68 


-0.06(0.04) 


-0.05(0.08) 




Mg/F 


G obs 


+0.44(0.09) 


+0.83(0.14) 






Mg/F 


e cor (EOO) 


1 n on/n i i \ 
+0.29(0.11 j 


+0.19(0.14J 






Mg/F 


ecor (Ell) 


+U.2y(U.ii ) 


+U.2D(U.iDj 






Mg/F 


e]cor (SOO) 


1 n 1 n/n i i \ 
+0.19(0.11 j 


+0.09(0.12 ) 






Mg/F 


e]cor (Sll) 


1 n 1 n/n i i \ 
+0.19(0.11 ) 


+0.20(0.1 ( ) 




[Ti/Fe 


obs 


< -0.09 


< -0.45 


+0.70(0.09) 




S/Zn 


obs 


+0.09(0.07) 


+0.01(0.12) 


+0.12(0.14) 




S/Fe]< 


ibs 


+0.34(0.07) 


+0.76(0.12) 


+0.62(0.15) 




S/Fe]< 


or (EOO) 


+0.19(0.09) 


-0.02(0.16) 


+0.25(0.19) 




S/Fe]< 


or (Ell) 


+0.19(0.09) 


-0.10(0.17) 


+0.25 0.19) 




S/Fe< 


or (SOO) 


+0.08(0.09) 


-0.30(0.17) 


+0.11(0.19) 




S/Fe]< 


or (Sll) 


+0.08(0.09) 


-0.30(0.17) 


+0.11(0.19) 


[N/S]obB 


-0.83(0.14) 


< -0.63 


< -0.85 




N/Si< 


abs 




< -0.48 


< -0.64 




N/Sil 


,or (EOO) 




< -0.66 


< -0.66 




N/Si] 


,or (Ell) 




< -0.72 


< -0.66 




N/Si 


,or SOO) 




< -0.90 


< -0.70 




N/Si] 


,or (Sll) 




< -0.87 


< -0.70 




Cl/S], 


)bs 




> -0.11 





Galactic stars the a-elements track more or less each other 
(within ±0.10 dex). Thus, by comparing a-elements of 
different dust depletion levels, their relative ratios when 
differing from the solar value are indicative of the pres- 
ence of dust. In the DLA toward QOlOO+13 we measure 
[S/Mg] = -0.10 ± 0.12, where Mg is relatively strongly 
deplet ed in th e Galactic ISM and S is a non-refractory ele- 
ment l)Savage fc SembachllQQfil) . This ratio is in agreement 
with the measured [Zn/Fe] ratio, both of them suggest a 
small amount of dust in this DLA system. In the DLA to- 
ward Q1331+17 we measure [S/Si] = +0.16 ±0.12. Since 
Si is only a mildly refractory element ijSavaee fc SembachI 
the observed S overabundance relative to Si indi- 
cates important amount of dust. We also have the mea- 
surement of [S/Mg] = -0.07 ± 0.18 in this DLA sys- 
tem. This abundance ratio rather indicates no presence 
of dust, but we ascribe this result to the weak accuracy 
of the Mgli column density measurement, which may be 
slightly overestimated due to some blends in the Lya for- 
est (see Paper I and Sect. l3.2|l . Finally, in the DLA toward 
Q2231-00 we measure [S/Si] = +0.21 ± 0.16 and [Si/Ti] 
= -0.25 ± 0.09. The [S/Si] ratio is in agreement with the 



conclusions derived from the [Zn/Fe] ratio, which high- 
lights the presence of dust in this DLA system. On the 
other hand, the measured undersolar [Si/Ti] ratio, while 
Ti is a strongly refractory element, does not confirm these 
findings. At a first glance, the undersolar [Si/Ti] ratio can 
only be the result of an overestimation of the Ti ii column 
density. 

5.2. Dust corrections 

Three main approaches can be used to circumvent the 
problem of dust depletion in the studies of DLA abun- 
dance patterns. The first one consists in considering 
only the DLA s ystems with wo or low dust dep l etion 
[Zn/Crl < 0.3 llPettini et a,].l l200fll: iMolaro et all 120001: 
ILodcz et al.1 l2002|) . Indeed, if the amount of dust is 
small, the dust depletion levels are negligible. This is 
valid for the DLA toward QOlOO+13, which has a low 
[Zn/Cr] +0.12 ± 0.02 ratio. However, for the two 
other DLAs, other solutions have to be found. The sec- 
ond approach consists in trying to quantify the dust de- 
pletion effects and provide dust corrections to the mea- 
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sured abundances. Several authors have considered dif- 
ferent met hods to compute the dust deplet ion correc- 
tions (e . g. IVladilol Il99a ISavadio et al.ll200(]|) . Recently, 
IVladild |2002allJ) elaborated a new and more complete 
method for correcting for dust depletion. Finally, the 
third approach consists in focusing on non-refractory and 
mildly refractory elements, such as N, O, S, Ar, and 
Zn. In this way one has directly access to the intrinsic 
abuiidances of DL As ('Centurion et al.' !200fltlMolaro et al.l 
l200nHDessauges-Zavadskv et aL.2(mii} . This approach to 
tackle the problem of dust depletion is by far the most ac- 
curate one. It is applicable to all DLAs whatever their 
amount of dust, and is independent of any assumption on 
the properties of dust in DLAs and of errors that may be 
introduced by a complex dust correction procedure. 

In our study we considered the second and the third 
approach to circumvent the problem of dust depletion. In 
Table [S] we summarize all the absolute abundances, and 
in Table |^ we summarize the relative abundances which 
are further used in the paper. Both the observed and 
dust-corrected values are reported. The dust cor rections 
were performed using the method developed by IVladild 
(,2002a^b '). This method groups together several dust cor- 
rection models based on different assumptions. They are 
labeled EOO, Ell when one assumes that the intrinsic 
[Zn/Fe] ratio is equal to -1-0.10 dex and SOO, Sll when the 
intrinsi c [Zn/Fe] ratio is equal to -f 0.00 dex (see more de- 
tails in lVladilcll2002a|) . The absolute abundances, [X/H], 
were obtained by summing the contributions of all the 
components reported in Tables ISHEl The relative abun- 
dances, [X/Y], were computed by considering only the col- 
umn densities of the components detected both in the X 
and Y profiles. In this way one avoids an overestimation of 
abundances derived from strong metal-line profiles relative 
to the abundances derived from weaker metal-line profiles 
in which generally only the strongest components are de- 
tected. In the case of very weak lines, like the Ti ii lines for 
instance, one can indeed underestimate the [X/Fe] ratios 
by up to 0.3 — 0.4 dex by considering the total Fe abun- 
dance generally derived from strong Fe ii lines. In the three 
DLAs studied this effect is particularly important in the 
DLA toward Q2231— 00 which show complex metal-line 
profiles, with a large number of components. The same 
approach has already been used in Paper I. 

5.3. Intrinsic abundance patterns 

In this Section we describe the derived intrinsic abun- 
dance patterns of the DLAs toward Q0100-hl3, Q1331+17 
and Q2231— 00. Their interpretation will be presented in 
Sect. El To present the abundance patt erns we use the 
same t ype of diagrams as the ones used bv lProchaska et all 
l)2003|) . The dust-corrected elemental abundances are plot- 
ted on a logarithmic scale, where hydrogen is defined to 
have e(II) = 12 and, in a general way, for any element 
X, e(X) = log(X/II) -|- 12, a nd are compared with the 
solar abundance pattern from ICrevesse fc Sauval (,1998i) 



(shown by the solid line) scaled to match the observed 
sulfur metallicity of each DLA system. This way of show- 
ing the elemental abundances as a function of the atomic 
number, Z, has the advantage that we see in the same plot 
all the measured abundances in a system, and this allows 
us to directly identify possible deviations from solar val- 
ues. 

Figures El EH and ^2 show the nucleosynthetic pat- 
terns of the DLAs toward Q0100-fl3, Q1331-hl7 and 
Q2231— 00, respectively. At a first glance, the abundance 
patterns of these high redshift galaxies resemble that of 
the solar neighborhood indicating that their nucleosyn- 
thetic enrichment histories are not too dissimilar from our 
Galaxy. However, at closer inspection, one notes several 
important differences. 

5.3.1. Q0100-M3, zabs = 2.309 

The DLA toward Q0100-hl3 (Fig. shows a shght en- 
hancement of the a-elements S and Mg abundances rela- 
tive to the Zn and Fe-peak element abundances, namely 
[S/Zn] = -h0.19±0.09 and [Mg/Fe]cor = 4-0.29±0.11. This 
a/Fe-peak enhancement is suggestive of an enrichment by 
massive stars. Indeed, the a-elements are produced in less 
than 2 x lO"^ yrs by Type II supernovae (SNe) resulting 
from massive stars and the Fe-peak elements are mainly 
produced by Type la SNe on longer timescales, between 
10^ - 10^ yrs. 

The dust-corrected ratio of the two a-elements Mg and 
S, [Mg/S]cor — -fO.lOibO. 12, shows that these two elements 
closely track each other in the DLA. On the other hand, 
the a-element Ar is slightly underabundant relative to S, 
as noted in Sect. ^ The Fe-peak elements Fe, Ni and Cr 
show solar values one relative to the other, as observed in 
the Galactic stars with similar metallicities (see Table (HJ. 

In this system, one also observes the odd-even effect, 
which corresponds to an underabundance of odd-Z ele- 
ments relative to the even-Z elements of the same nucle- 
osynthetic origin. Indeed, we obtained [P/S] = — 0.37 ± 
0.12 at [P/H] = —1.85 ± 0.13, which even shows evidence 
for an enhanced odd-even effect. This value is similar to 
the [P/Si] = -0.21 at [P/H] = -1.16, [P/Si] = -0.40 
at [P/H] = -2.30 and [P/Si] = -0.30 a t [P/H] = -1.20 
values measured toward other DL As bv iLevsh akov e t all 
l|2002l) . iMolaro ct al. (200l]) and lOutram et al. [ IdM), 
respectively, and is indicative of massive supernovae. 
Phosphorus abundance measurements in DLAs are very 
important to understand the nucleosynthesis of this el- 
ement, because no P abundance measurement exists 
in the Galactic stars. The observed odd-even effect in 
the P/S ratio is in agre ement with the findings of 
iGoswami fc PrantzosI l|200n>) . 

Finally, the DLA shows undersolar [N/S] ratio of 
-0.83 ± 0.14 at [S/H] = -1.48 ± 0.11. N is an element 
of particular interest, since it has a complex nucleosyn- 
thetic origin. Synthesized in the CNO cycle of stars, N 
is either a secondary element in the sense that it is pro- 
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Fig. 9. The nucleosynthetic abundance pattern of the DLA at Zabs = 2.309 toward QOlOO+13. The dust-corrected 
elemental abundances are represented on a logarithmic scale, where hydrogen is defined to have e(H) = 12 and in a 
general way, for any eler nent X, e(X) = log(X/H) + 12. They are compared with the solar abundance pattern from 
ICrevesse fc Sauyal l|l998() shown by the sohd line and scaled to match the obse rved sulfur metallicity, [S/H] = —1.44, of 
the DLA system. The dust corrections were obtained via the lVladilol l)2002albf) method. The filled circles correspond to 
the dust-corrected abundances obtained with the EOO dust correction model and the open squares to the dust-corrected 
abundances obtained with the Ell dust correction model. 



duced in proportion to C and O originally present in the 
star, or a primary element in the sense that it can be pro- 
duced starting from C and O manufactured by the star "in 
situ". In the N/a versus a/H diagram, the secondary N 
over a-elcmcnt ratio is expected to increase steeply with 
the increasing metallicity, while the primary N over a- 
element ratio remains constant when the metallicity in- 
creases (jTalbot fc Arnettlll974|) . The progenitors of pri- 
mary N have not bee n fully constrained yet. Th e most re- 
cent stellar models of lMevnet fc MaedeJ l)2002|) including 
the stellar rotation indicate that the intermediate-mass 
stars are the main producers of N. The [N/S] ratio reached 
in the DLA toward QOlOO-l-13 is a "high" value, in the 
sense that it is close to the primary N "plateau" at [N/S] 
~ —0.75. It is at the upper end of the range of N/S values 
measured in DLAs and in good agreement with the typi- 
cal values observed in Hll regions at similar me tallicities 
(e.g. Ilzotov fc Thuanlligggt IPilvuein et alll2002() . 

5.3.2. Q1331-hl7, z^bs = 1-776 

The DLA toward Q1331+17 (Fig. ^ has almost solar 
abundance pattern. Very few differences are observed. 
Indeed, the Fe-peak elements Fe, Ni and Cr closely track 
each other, the two a-elements S and Si also (see TablelH)), 
and the a-element over Fe-peak element ratios — Si/Fe, 
S/Zn and S/Fe — are all solar. Only the Mg/Fe ratio is 



slightly oversolar (see Table EJ. We would like to note 
that the SOD and Sll dust correction models of I Vladiiol 
l|2002al) which assume a solar ratio of Zn/Fe, lead to 
[S/Fe]cor = -0.30 ±0.17 (see Table EJ. Such an under- 
solar a-element over Fe-peak element ratio is neither ob- 
served in Galactic stars nor predicted by chemical evolu- 
tion models. Thus, this provides evidence that [Zn/Fe] > 
(intrinsically) in this DLA system. Ti exhibits a large un- 
derabundance relative to the Fe-peak elements and this, 
although Ti is an a-element, and even after having ap- 
plied the dust corrections. At a first glance, the only way 
to explain this observation is that we have underestimated 
the Ti upper limit. 

We measure the abundances of P and Mn, two odd-Z 
elements. The odd-even effect is less marked in this system 
than in the DLA toward Q0100-hl3. We get [P/Si]cor = 
—0.20 ± 0.11 when the EOO dust correction is applied 
and [P/Si]cor = -0.26 ± 0.12 when the Ell dust cor- 
rection is applied at [P/H] = -1.42 ± 0.13, and [Mn/Fe] 
= -0.15±0.04 (without dust corrections, and -0.07±0.05, 
— 0.37± 0.05 when the EOO, Ell dust corrections are ap- 
plied, respectively) at [Fe/H]cor = —1.23 ±0.12. 

We have detected the CI i A1347 fine in this DLA sys- 
tem, which allows us to measure the Cl*^ column den- 
sity and to provide a lower limit on the total CI abso- 
lute abundance. CI is also an odd-Z element. This is the 
third CI abundance upper limit derived in a DLA sys- 
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Fig. 10. The nucleosynthetic abundance pattern of the DLA at Zabs = 1.776 toward Q1331+17. The solar values were 
scaled to match the observed sulfur metallicity, [S/H] = —1.26. Same as Fig. El 



tem ^ the other two were ob tained bv lLedoux et all l|2002(l 
and lProchaska et all l)2003(l . CI abundance measurements 
in DLAs are very important to constrain the CI nucle- 
osynthesis in the low metallicity regime, in particular be- 
cause no CI abundance measurement exists in Galactic 
stars. CI measurements are available in planetary neb- 
ulae of the Galact ic disk (PNe; e.g. ICosta et alJ Il996t 
Kwitter et all l2003l) and in the Galactic disk ISM (e.g. 
Weltv et alJll999|) . Oversolar [Cl/S] ratios between +0.09 
and -1-0.5 dex are observed in Type II PNe, and under- 
solar [Cl/S] ratios of —0.6 and —0.2 dex are measured in 
the cold and warm Galactic disk, respetively, suggesting 
that CI is dust depleted. In the DLA studied here we de- 
rive a lower limit on the [Cl/S] ratio of > -0.11 at [Cl/H] 
> - 1.37, and in the D LA studied bv lLedoux et all (2002^ 
and lProchaska et alJ (^.2003) [Cl/S] > -0.33 and > -0.68, 
respectively. These three lower limits measured in DLAs 
suggest that the Cl/S ratios likely are solar to oversolar 
similarly to what is observed in PNe, thus leaving very lit- 
tle space for an o dd-even effect for CI. The ch emical evo- 
lution models of Goswami & Prantzos ('2000') computed 
with the Woosle v fc Weaver. (J,2£5) yields predict a solar 
Cl/S ratio, and hence a negligible odd-even effect for CI. 

IProchaska fc Wol^ l|l999ll also reported the detection 
of Cii* in the DLA toward Q1331-t-17. The presence of 
both strong C II* and CI i absorptions suggests that the 
gas resides in a cold neutral medium, characteristic of 
highly depleted gas in the Milky Way. This DLA precisely 
shows an important dust depletion level, with [Zn/Fe] 
= +0.75 ± 0.05. Furthermore, the observation of signifi- 
cant CI I requires at least a modest molecular hydrogen 
fraction. Indeed, in regions where II2 is optically thick. 



CI II reacts rapidly with II2 to form HC1+, which in turn 
leads to CI i and H i. Because the conversion of CI 11 to CI i 
is faster than the photoionization of CI i, CI is primarily 
neutral in regions where H-^ is abund ant but otherwise is 
primarily ionized ij.Tura fc Yorklll978l) . 

Finally, the derived upper limits on [N/S] and [N/Si] 
are at the upper end of the range of N/a values measured 
in DLAs (see Table . They are close to the primary N 
"plateau" . 

5.3.3. Q2231-00, Zabs = 2.066 

The DLA toward Q2231-00 (Fig.jnj shows a solar [Si/Fe] 
ratio, but enhanced [S/Zn,Fe] and [Ti/Zn,Fe] ratios (see 
Table O . The S abundance was derived from a single S 11 
line located in the Lya forest (see Sect. 13. 3t . and we believe 
that it may be partly blended with Hi clouds. The com- 
parison of the dust-corrected absolute abundances of Si 
and S, two a-elements, shows further evidence for a slight 
overestimation of the abundance of S, since [S/Si]cor — 
-1-0.19 ± 0.16 when the EOO dust correction is applied. 
Although the adopted large error on the S abundance mea- 
surement takes into account this slight overestimation, in 
Fig. ^3 we scale the solar values to match the dust cor- 
rected silicon metallicity, instead of the observed sulfur 
metallicity as done in Figs. 1^ and IIUI The Ti abundance 
is also very difficult to measure, thus it is very likely that 
it is overestimated, as already suspected in Sect. 15.11 

The Fe-peak elements, Fe, Ni and Cr, track very well 
each other, as observed in the Galactic stars with similar 
metallicities (see Table [HI. With the abundance measure- 
ment of the Fe-peak element Mn, we highlight a rather 
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pronounced odd-even effect in this DLA system, [Mn/Fe] 
= -0.16 ± 0.05 at [Fe/H]cor = -0.83 ± 0.14. Another 
odd-Z element detected in this system is P. Only an 
upper limit on the [P/Si] abundance ratio is derived, 
[P/Si]cor < +0.34, which is potentially in line with an 
odd-even effect for P. 

The obtained [N/Si,S] upper limits (see Table El are 
again "high", close to the primary N "plateau". We do 
not believe that the observed N i lines are heavily blended 
with Lya forest lines. 

6. Determination of the star formation histories 

Until now the DLA galaxy population has been ana- 
lyzed as a whole. Several chemical evolution models were 
constructed in order to interpret the abundance pat- 
terns observed in DLAs as an ensemble, considering them 
as an evolutionary sequence, i.e. obje cts caught in a 
differen t phase of their evolution (e.g. 'Matteucci et al.| 
1997: .Timenez. Rowen Ma,tteuccii1 99fl:..Hou et al. 2001i 
MathHn et alJl200ll) . However, several pieces of evidence 
— like the low redshift deep imaging revealing a variety 
of m orphological types be l onging to the DLA population 
(e.g. iLe Brun et all Il997t iNestor et aLll2002|) . the large 
scatter in the a over Fe-peak element abundance ratios at 
a given metallicity and the large scatter observed in the 
metallicities — indicate that the DLAs trace galaxies with 
different star formation histories. Some may have formed 
stars on timescales similar to that of the early Milky Way, 
while others apparently did so more slowly or intermit- 
tently, so that the Fe-peak elements could catch up with 
the a-elements. 



Thanks to the large number of elements detected in 
the DLA systems studied in this paper, we are for the 
first time in a very appropriate situation to analyze these 
high redshift galaxies individually. The different elemental 
abundances at disposal allow us, in particular, to constrain 
the star formation history (hereafter SFH, and SF for the 
star formation) and the age of the associated galaxies by 
means of a direct comparison with a grid of chemical evo- 
lution mo dels for spiral and dwarf irregular galaxies (see 
ICalura eTal. 20031 rather than just by means of a com- 
parison with the solar abundance pattern as presented 
in Sect. 15.31 We present below the approach which has 
been applied. The derived results for the DLAs toward 
QOlOO+13, Q1331+17 and Q2231-00 are summarized in 
Table Cni 

6.1. Chemical evolution models 

A chemical evolution model allows one to follow in detail 
the evolution of abundances of several chemical species, 
starting from the matter reprocessed by the stars and 
restored into the ISM through stellar winds and super- 
nova explosions. Here we briefly summarize the major in- 
gredients of the chemical evolutions models used in this 
work which we identify as "spiral" and "dwarf irregular" 
models according to the type of galaxies they do match 
best. These raodels are the same as the models used by 
'Calu ra etall l)2003|) in their study of the DLA nature. A 
de tailed description of the "sp iral" model can be found 
in lChiappini et all lll997Ll200ll) and of the " dwarf irregu- 
lar" inodel i n iBradamante et alJ l|l998') and^ R,ecchi et alJ 
l|2nfllLl2002ll . In both models no instantaneous recycling 
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DLA system 




Model 




Star formation characteristics 
V [Gyr-i] / At [Gyr] / % [Gyr] 






Age 

[Gyr] 


QOlOO+13, 2abs 


= 2.309 


Spiral 




solar neighborhood: R = S kpc 




2.5 


0.25 ±0.20 






Dwarf irreg 


ular 


bursting SF: 0.9 / 0.07 / 2.0 


(2.35 


~ 2.47) 


0.05 - 0.20 


Q1331 + 17, 2abs 


= 1.776 


Spiral 




outer regions: _R = 16 kpc 




2.8 


1.5 ±0.4 






Dwarf irreg 


ular 


continuous SF: 0.03 / 13 / ... 




- 10 


> 3.5 


Q2231-00, 2abs 


= 2.066 


Dwarf irreg 


ular 


bursting SF: 4.2 / 0.10 / 0.5 


(2.10 


-2.16) 


0.05-0.15 



Zf in brackets — Values obtained under the assumption that the galaxy has undergone a single burst of SF 



approximation is adopted, namely the stellar lifetimes are 
properly taken into account. 

6.1.1. "Spiral" model 

The spiral galaxies are a ssumed to form as a result of 
two main infall episodes (|ChiaDDini et alJll997l) . During 
the first episode the halo forms and the gas shed by the 
halo rapidly gathers in the center leading to the forma- 
tion of the bulge. During the second episode, a slower 
infall of external gas gives rise to the disk with the gas 
accumulating faster in the inner than in the outer re - 
gion ("inside-out" scenario. iMatteucci fc FrancoisHlQS^ . 
The process of disk formation is much longer than the 
halo and bulge formation, with timescales varying from 
~ 2 Gyr in the inner disk to 7 Gyr in the solar re- 
gion and up to 13 Gyr in the outer disk. In particular, a 
timescale of 7 Gyr in the solar neighborho od is required to 
fit the G-dwarf metallicity distr ibution ( Chiappini et alJ 
Il997t iBoissier fc Prantzojll999(l . The adopted SFR ex- 
pression is: 



|2(fc-l) 



.a{R,t 



Gal) 



C7{R,t) 



-'at,AR,t) (1) 



where v is the star formation efficiency, namely the in- 
verse of the typical timescale for star formation, a{R,t) 
is the total (gas + stars) mass surface density at a radius 
R and time t, a{RQ,t) is the total mass surface density 
in the solar region, and <Jgas{R,t) is the gas mass sur- 
face density. For the gas density exponent k a v alue of 
1.5 has been assumed by IChiappini eFal] l)l997|) m or- 
der to ensure a good fit to the observational constraints 
in the solar vicinity. The star formation efficiency is set 
to v ~ 1 Gyr~^, and becomes zero when the gas sur- 
face density drops below a certain critical threshold. A 
threshold densit y ath ~ 7 Mp^ PC~^ is a dopted in the disk 
as suggested bv lChiappini et aP ill 9971). T he initial mass 
function (IMF) is taken from'Scalo' ('1986') The galactic 
wind is considered not efficient (:Matteucciii200lh . All the 
input parameters of this "spiral" model were constrained 
in order to successfully reproduce the observed properties 
first in the solar vicinity, such as the solar abundances, 
the G-dwarf metallicity distribution (see above) and the 
abundance ratios, and then at the global scale, such as the 
radial gas and stellar density profiles and the abundance 



gradients of various elements observed in the stars of the 
Milky Way 

Due to the "inside-out" scenario and the threshold 
density adopted in the "spiral" model, the SFH is differ- 
ent at different galactocentric radii. Indeed, in the inner 
regions of the disk where the rate of accretion of mat- 
ter onto the disk is fast, a high gas surface density, well 
above the critical threshold, is reached and maintained 
during a long period, and hence the SF is almost continu- 
ous through the galaxy lifetime. On the other hand, in the 
outer regions of the disk (radius R > 12 kpc) where the 
rate of accretion of matter onto the disk is slow, the SFH 
proceeds in a gasping way, due to the fact that in these re- 
gions the gas is always close to the critical threshold. As a 
consequence, the external regions of the disk look like the 
Magel lanic irregulars w here this kind of SF seems to take 
place ((Tosi et al.ll99lj) . These different star formation his- 
tories at different galactocentric radii are associated with 
different abundance patterns which can be compared with 
observations. 

When comparing the "spiral" model with our obser- 
vations of DLA systems, we run a set of models each of 
them corresponding to a different galactocentric radius 
R. This radius represents the position at which the QSO 
line of sight crosses the disk of the observed DLA galaxy. 
Finally, we need to constrain the formation redshift, Zf, 
of the "spiral" model, which is the single free parameter. 

6.1.2. "Dwarf irregular" model 



We ad opt a model based on the work of lBradamante et alJ 
lll998l) . in which the dwarf irregular galaxies form ow- 
ing to a continuous infall of pristine gas with an infall 
timescale of 0.5 Gyr, until a mass of 10^ M0 is accumu- 
lated. The evolution of dwarf irregular galaxies is assumed 
to be characterized by a bursting star formation history. 
This particular model was built to reproduce the chemi- 
cal properties of the local blue compact galaxies (BCGs). 
The parameters which need to be defined in this model 
are the number of bursts of SF that the galaxy has un- 
dergone, and for each burst, the star formation efficiency, 
v, (same definition as in the "spiral" model), the burst 
duration. At, and the time of occurrence of the burst, th, 
i.e. the time necessary for the infall of pristine gas before 
the SF starts. Hence, the star formation rate in the "dwarf 
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irregular" model can proceed either in short bursts of a 
duration from 10 to 200 Myr separated by long quiescent 
periods, or at a low regime but continuously, namely in 
one or two long episodes of inefficient SF lasting between 
3 and 13 Gyr. In this model, the adopted SFR is: 

^(t) ^ lyagasit) (2) 

where agas{t) represents the gas density at time t and v 
is the star formation efficiency. 

The dwarf irregular galaxies are particularly sensi- 
tive to outflows resulting from the energy injection from 
the SN explosions. Indeed, galactic winds develop when 
the thermal energy of the gas equals its binding energy 
ijMatteucci fc Tornamblll987|) . In the "dwarf irregular" 
model used in this paper, the rate of gas loss via galac- 
tic winds is assumed to be proporti onal to the SF, an d 
the recent chemodynamical results of lRecchi et al ](|200ll), 
suggesting that the wind is differential, are adopted. This 
implies that the ejecta of Type la SNe and intermediate- 
mass stars are lost from the parent galaxy more easily 
than Type II SN ejecta, so that the Fe and N ejection 
efficiencies are larger than the ejection efficiencies of a- 
elements (e.g. O, Mg). This is due to the fact that the 
Type la SNe can transfer into the ISM more energy than 
the Type II SNe, since they explode in a n already r arefie d 
and hot medium. The IMF is taken from lSalpeteJ l)l955|) . 

For the comparison of the "dwarf irregular" model with 
our observations of DLA systems, we assume that these 
high redshift galaxies have undergone a single burst of SF. 
This can be justified by the fact that the DLA galaxies 
are probably young galaxies. We thus need to essentially 
constrain three free parameters^ the burst star formation 
efficiency, the burst duration. At, and the time of oc- 
currence of the burst, t^. 

6.1.3. Stellar yields 

We describe here the nucleosynthesis prescriptions 
adopted in both the "spiral" and " dwarf irregu- 
lar" in odels. They include the yields of iNomoto et alJ 
lll997al) for massive stars (M > lOM©), the yields 
of Ivan den Hoeck fc Groenweeenl l)l997j) for low- and 
in termediate-ma, s s stars (0.8 < M/Mq < 8) and the yields 
of lNomoto et all l)l997 b') for Type la SNe (model W7). For 
Zn and Ni we consider the following specific prescriptions. 

The nucleosynthesis of Zn is a debated issue, since 
what process is making Zn is rather uncertain. According 
to stellar models, the production of Zn can ensue via s- 
processes in low- and high-mass stars during He-core burn- 
ing as well as during explosive nucleosynthesis events oc- 
curring in Type la and II SNe t,Matteucci et al. .1993. and 
references therein). This uncertainty in the Zn production 
has repercussions on the reliability of Zn yields. In this pa- 
per we use the Zn yields adopted bv lCalura et al. ( 2003) , 
which were estimated by extrapolating the Galactic star 
abundance observations in the solar neighborhood to the 
high redshift Universe, assuming Fe and Zn track each 
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Fig. 12. Observed and predicted [Ni/Fe] versus [Fe/H] 
distributions in the solar neighborhood for two different 
nucleosynthetic prescriptions for the production of Ni. The 
dashed line corresponds to the model with /? = 0.0065 
and the solid line to the model with (3 — 0.0049, where 
/9 is a multiplicative factor between the Fe a nd Ni yields 

isee text). The solid squares are the d ata fromlRvan et alJ 
1996), the open circles the ones fromlChen et alJ ll2000l) . 
the open triangles the ones from Jehin ct al." (1999), the 
filled triangles the ones from iSteohens, (■1999i') . the filled 
circles the ones from lMcWillian et al .Illl995l). and the open 
squares the ones from lGratton & SnedenI l|l99lh . 



oth er as observed in the metallicity range where t he DLAs 
lie l|Prochaska et alJl200(ll: iMishenina et al.ll2002l) . 

The Ni yields have also some difficulties to reproduce 
the Galactic star abundance observations. Therefore, sim- 
il arly to Zn, we assumed for the Ni production the results 
of lMatteucci et all l)l993) . in which a good fit to the solar 
abundance of Ni is found if the bulk of its production is 
ascribed to Type la SNe. However, since a non-negligible 
fractio n of Fe also comes from Type II SNe l|Nomoto et alJ 
Il997al) and since the abundance of Fe seems to vary in 
lockstep with that of Ni in the solar vicinity, the Type 
II SNe should also produce some Ni. In our models the 
amount of Ni produced through explosive nucleosynthe- 
sis in massive stars scales with the Fe yields, according to 
^Ni = P X Ypc, where /? represents a multiplicative fac- 
tor. We have run a chemical evolution model for the solar 
vicinity varying /3 in order to reproduce the [Ni/Fe] versus 
[Fe/H] distribution observed in Galactic stars of different 
metallicities. Figure IT^ shows this [Ni/Fe] versus [Fe/H] 
distribution compared with the predictions for the solar 
neighborhood model, when two values of P are adopted 
(dashed line: /3 = 0.0065; sofid line: /3 = 0.0049). The 
choice of /3 = 0.0049 gives a satisfactory agreement be- 
tween the predictions and the observations, in particular 
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in the metallicity range occupied by the DLAs ([Fe/H] 
< —2). For Ni produced in Type la SNe, we assume a 
constant value Mm; ^ 2.3 x 10"^ Mfr,. 



6.2. Comparing the DLA abundances with the models 

To be able to constrain the star formation history of DLAs 
from the comparison of their abundance patterns with 
chemical evolution models, we first need to understand 
how the abundance patterns can provide information on 
the star formation history of galaxies. In few words, the 
absolute abundances depend on the model assumptions 
(i.e. the SFH), whereas the relative abundances depend 
only on the nucleosynthesis, the stellar lifetimes and the 
IMF. Relative abundances can therefore be used as cos- 
mic clocks if they involve two elements formed on differ- 
ent timescales, like it is the case of the a-element over 
Fe-peak element ratios and the N over a-element ratios. 
As a consequence, these abundance ratios when exam- 
ined together with the absolute abundances [Fe/H], or any 
other metallicity tracer such as [Zn/H], and [a/H], respec- 
tively, completely determine the .star formation history of 
a galaxy. On the other hand, when these abundance ratios 
are examined as a function of the redshift, they provide 
constraints on the age of a galaxy, defined as the epoch 
at wh ich the galaxy has started to form stars I Matteucc3 

The relative abundances that we have at disposal in 
the DLA systems studied here are [Si/Fe], [S/Zn], [S/Fe], 
[Mg/Fe], [Ni/Fe], [N/Si] and [N/S]. For each DLA individ- 
ually, we consider all of them at once to first constrain the 
SFH and then the age of the associated galaxy by means 
of a detailed comparison with a grid of chemical evolution 
models. This grid includes the "spiral" model computed at 
different galactocentric radii from i? = 2 to 18 kpc and the 
"dwarf irregular" model computed for a single burst of SF 
with different characteristics defined by varying the three 
free parameters v, At and ti,. The adopted cosmology is 
Ho = 65 km s~i Mpc'^ (and hence h = 0.65), J7m = 0.3 
and r^A — 0.7. 

To determine the best chemical evolution model repro- 
ducing the abundance patterns observed in each DLA sys- 
tem by taking into account all the different informations 
available on the system, we use a statistical test when the 
best solution cannot be clearly identified by eye. This test 
consists first in determining the minimal distance between 
the data point and the curve of a given chemical model 
in each abundance diagram at disposal. For this purpose, 
we consider the 1 a error on the data point (or more pre- 
cisely the covariance matrix of the 1 a measurement er- 
ror), and we derive this minimal distance by computing 
the distances di from the data point to the points defin- 
ing the considered model curve and by looking for the 
di for which the di/ai ratio is minimal. Second, once the 
minimal distances for all the abundance diagrams consid- 
ered in each system and for a given model are derived, we 
compute their weighted mean. Finally, the comparison of 



the weighted means obtained for different chemical mod- 
els determines the best chemical evolution model which 
represents the data points and thus the DLA galaxy. The 
upper and lower limits are not taken into account in this 
test. 



6.2.1. Q0100-M3, Zabs = 2.309 

To determine the nature of this DLA galaxy, we have five 
different abundance ratios at disposal. We can thus try to 
constrain the best chemical evolution model which repro- 
duces this DLA galaxy by taking into account the informa- 
tion provided by the following five diagrams: [S/Zn] ver- 
sus [Zn/H], [S/Fe] versus [Fe/H], [Mg/Fe] versus [Fe/H], 
[Ni/Fe] versus [Fe/H] and [N/S] versus [S/H]. For the ab- 
solute and relative abundances of refractory elements we 
consider the dust-corr ected va, l ues as der ived with the dust 
correction method of IVladiiol l)2002albl . and see Tables |H1 
and ISJ. In this DLA system, however, the dust correc- 
tions are small, since the [Zn/Fe] ratio is close to solar 
(= +0.25 ±0.04). 

We considered the "spiral" model at the galactocen- 
tric radii i? = 6, 8, 12, and 18 kpc. The "spiral" model at 
R = 8 kpc best reproduces the five abundance ratios mea- 
sured in this high redshift galaxy (see the thick dashed- 
dotted curves in Fig. I13|) . We might be satisfied with this 
solution, since the model curves are in agreement with the 
data points within less than 1 ct in four out of the five di- 
agrams. In the case of the Ni/Fe ratio, the model is in 
agreement only within 2 — 3 ti. But, to have a complete 
picture and to check the uniqueness of the derived solu- 
tion, we investigated whether a "dwarf irregular" model 
can also correctly reproduce the data points. We consid- 
ered a "dwarf irregular" model with a single burst of SF 
and explored the following values for the three parame- 
ters characterizing the burst: = 0.5, 0.7, 0.9, 1 Gyr~^, 
At = 0.05, 0.07, 0.1, 0.2 Gyr and 4 1, 1.1, 1.5, 2 Gyr. 
The "dwarf irregular" model which best reproduces the 
data points has one burst with a star formation effi- 
ciency ly = 0.9 Gyr~^ and a short duration At = 0.07 
Gyr (see the thin solid curves in Fig. I13|l . The burst 
occurs after th ^ 2 Gyr of continuous infall of pristine 
gas. These "dwarf irregular" model characteristics very 
well match the starburst parameters dete rmined for the 
BCGs (e.g. iLanfranchi fc Matteuccil Eo03l) . In summary, 
two chemical evolution models yet tracing different star 
formation histories correctly reproduce the abundance ra- 
tio versus metallicity distributions observed in the DLA 
toward QOlOO+13. The "spiral" model at i? = 8 kpc has 
a continuous SFH, whereas the "dwarf irregular" model 
constrained by the observations has a bursting SFH (see 
Sect. I6.1|l . The identification of the best model among 
these two is very difficult, because at this intermediate 
metallicity both models yield very similar chemical abun- 
dance patterns. The distinction is more pronounced at 
[Fe/H] < -2 or at [Fe/H] > -1. 
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Fig. 13. Observed and predicted abundance ratios versus metallicity for the DLA at Zabs = 2.309 toward QOlOO+13. 
The thick dashed- dotted curves correspond to the "spiral" model at i? = 8 kpc. The thin solid curves correspond to the 
"dwarf irregular" model with a single burst of SF occurring at = 2 Gyr and with a star formation efficiency v = 0.9 
Gyr~^ and a duration At = 0.07 Gyr. In this and all the following figures, the open squares represent the measured 
abundance ratio measurements, an d the filled trian gles and the filled squares represent the abundance ratios corrected 
for dust depletion according to the 

|viadi3 (HqoIS)^ 

method using the models EOO and Ell, respectively. 



We carried out the investigation of the age of the sys- 
tem, namely the most likely redshift of formation, 2:/, of 
the system, for the two models. For the "spiral" model 
at i? = 8 kpc, we analyzed the [Zn/H], [Ni/Fe], [Mg/Fe], 
[S/Zn], [S/Fe] and [N/S] versus redshift diagrams. The es- 
timated redshift of formation is between Zf = 2.4 — 3 with 
the best solution being at Zf 2.5 which is in excellent 
agreement with all the data points (see the thick dashed- 
dotted curves in Fig. I14|l . The DLA system is observed 
at Zabs — 2.309, hence Zf ^ 2.5 corresponds to an age 
of the DLA galaxy of 0.25 ± 0.20 Gyr. In the case of the 
"dwarf irregular" model, the determination of the age is 
more complex, because this model has several free param- 
eters. The parameter having the lowest impact is the time 
of occurrence of the burst, if,, according to the adopted 
definition of the age (see Sect. l6.2jl . In addition, ti, is diffi- 
cult to determine, whereas the star formation efficiency, v, 
and the burst duration. At, can be robustly constrained 
by the observations. It is the burst duration which has the 



largest weight in the determination of the age of the sys- 
tem. We satisfactorily reproduce the data points (within 
1 cr) for At values between 0.05 and 0.20 Gyr. Hence, we 
feel confident in concluding that the possible age of the 
system is between 0.05 — 0.20 Gyr. Given the assumption 
that the galaxy has undergone a single burst of SF, the 
corresponding redshift of formation is between 2.35 — 2.47. 

Both models point to a young age for the DLA galaxy, 
lower than 250 Myr. This age is consistent with the a- 
element over Fe-peak element enhancement observed in 
this DLA system (see Sect. I5.3ll . which indicates that the 
enrichment of this system is dominated by massive star 
products and requires a minimal contribution from Type 
la SNe. In addition, we saw in Sect. l5.3l that the N/S ratio 
in this DLA system is relatively high, very close to the pri- 
mary N "plateau", [N/S] = -0.83 ± 0.14. This N/S ratio 
and the estimated age for the DLA galaxy provide im- 
portant constraints on the N stellar progenitors. Indeed, 
such a high N/S value can be reached within less than 250 
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Fig. 14. Observed and predicted abundance ratios versus redshift for the DLA at Zabs — 2.309 toward QOlOO+13. The 
thick dashed- dotted curves correspond to the "spiral" model at i? = 8 kpc and a redshift of formation Zf — 2.5. For 
the definition of symbols, see Fig. 1131 



Myr, only if N is produced by intermediate-mass stars with 
masses between 5 — 8 Mq , which h ave lifetimes between 30 
and 70 Myr, and by massive stars. IChiaDpini et all l|2003|) 
recently showed that the same chemical evolution models 
as those used in this paper, b ut computed with t he re- 
cent pubhshed stellar yields of iMevnet fc Maede^ l)2002|) 
which take into account the effects of rotation in the stel- 
lar evolution, still reproduce such high N/S values within 
a timescale lower than 250 Myr. 

6.2.2. Q1331-hl7, Zabs = 1-776 

We have collected a lot of information on this DLA sys- 
tem. The observed abundance ratios at our disposal are 
[S/Zn], [S/Fe], [Si/Fe], [Mg/Fe], [Ni/Fe] and the upper lim- 
its on [N/S] and [N/Si]. This DLA system has a very high 
[Zn/Fe] = 0.75 ±0.05 ratio and exhibits one of the largest 
dust depletion level of any DLA. In Fig.^jwe can see the 
high differences between the observed and dust-corrected 
values. 

First, we examined the "spiral" model at the galacto- 
centric radii R = 8, 12, 16 and 18 kpc. The best model 



reproducing the data points in the [S/Zn] versus [Zn/H], 
[S/Fe] versus [Fe/H], [Si/Fe] versus [Fe/H], [Mg/Fe] versus 
[Fe/H] and [Ni/Fe] versus [Fe/H] diagrams is the "spiral" 
model of outer regions of the disk, at i? = 16 kpc (see the 
thick dashed-dotted curves in Fig. I15() . The model curves 
are in excellent agreement with the data points within less 
than 1 (7 in the five diagrams. Secondly, we investigated 
the "dwarf irregular" model to check the uniqueness of 
the solution. We considered the "dwarf irregular" model 
with a single burst of SF and explored a large range of val- 
ues for the three parameters characterizing the burst. We 
found out that the "dwarf irregular" model with a burst- 
ing SFH cannot reproduce the observations, which show 
solar a-element over Fe-peak element ratios, because the 
bursting SFH predicts a-enhanced patterns. On the other 
hand, the "dwarf irregular" model with a continuous SFH 
characterized by a single burst of SF with a star forma- 
tion efficiency in the range v — 0.02 — 0.09 Gyr~^, with 
the best value being v — 0.03 Gyr^^, and a duration over 
the whole Hubble time satisfactorily reproduces the ob- 
servations (see the thin solid curves in Fig. I15|l . 
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Fig. 15. Observed and predicted abundance ratios versus metallicity for the DLA at Zabs = 1-776 toward Q1331+17. 
The thick dashed- dotted curves correspond to the "spiral" model at i? = 16 kpc and the thin solid curves correspond 
to the "dwarf irregular" model with a continuous SF at an efficiency v = 0.03 Gyr~^. For the definition of symbols, 
see Fig. ^| 



Again none of these two models yet having a different 
star formation history — the "spiral" model of outer re- 
gions of the disk has a bursting SFH (see Sect. l6.l!T|) . while 
the "dwarf irregular" model constrained by the observa- 
tions has an inefficient continuous SFH — can be easily 
discarded, because they yield similar chemical abundance 
patterns due to the fact that they both are character- 
ized by a weak SF. The Mg/Fe ratio could be used as a 
discriminant between the two models (see Fig. I15|l . but 
there are substa ntial uncertainties in the Mg yields (e.g. 
IChiapDini et~anfl99 9') . 



We investigated the age of the DLA galaxy for the 
two models. For this purpose, we considered the [Zn/H], 
[S/Zn], [S/Fe], [Si/Fe], [Mg/Fe] and [Ni/Fe] versus redshift 
diagrams. The range of redshifts of formation determined 
for the "spiral" model at i? = 16 kpc is zj = 2.3 — 3, with 
the best solution being at Zf ^ 2.8 (see the thick dashed- 
dotted curves in Fig. ll6|) . For the "dwarf irregular" model 
with an inefficient continuous SF rate, the estimated red- 
shift of formation is large, roughly around z/ ~ 10 (see the 
thin solid curves in Fig. I16|l . The DLA system is observed 
at Zabs = 1.776, hence Zf ^ 2.8 and 2/ 10 correspond 
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Fig. 16. Observed and predicted abundance ratios versus redshift for the DLA at Zabs = 1-776 toward Q1331+17. The 
thick dashed- dotted curves correspond to the "spiral" model at i? = 16 kpc and a redshift of formation Zf — 2.8. The 
thin solid curves correspond to the "dwarf irregular" model with a continuous SF at an efficiency v = 0.03 Gyr^^ and 
a redshift of formation Zf — 10. For the definition of symbols, see Fig. 1131 



to an age of the DLA galaxy of 1.5 ± 0.4 Gyr and > 3.5 
Gyr, respectively. 

According to the recent WMAP results, the reion- 
ization seems to have occurred at a redshift consider- 
ably higher than what was thou ght previously, namely 
at Zrcion = 20+g° ijBennett et al.l [2003). Therefore, if the 
intergalactic medium has been reionized by the first stars 
formed in very young galaxies, a redshift of the order of 10 
could be reasonable for the appearance of the first galactic 
structures in the Universe. The derived solution that the 
associated galaxy with the DLA system at Zabs — 1.776 



toward Q133H-17 may be a dwarf irregular galaxy with 
a continuous SFH formed at Zf ~ 10 thus is possible. 
The solution that the associated galaxy is a spiral galaxy 
with the QSO line of sight crossing the outer regions of its 
disk at i? = 16 kpc and an age of 1.5 Gyr is also plausible. 
Both models point to a large age for this DLA galaxy. This 
refiects the long timescale necessary to accrete the gas as 
well as to reach the observed solar a-element over Fe-peak 
element ratios (see Sect. I5.3|l which require a substantial 
contribution from Type la SNe, releasing the Fe-peak el- 
ements only after 10* — 10^ yrs. But, both models fail in 



30 



M. Dessauges-Zavadsky et al.: A comprehensive set of elemental abundances in DLAs 



0) 



0.4 r 
0,2 - 


-0.2 r 
-0.4 q ^ 
-2.5 



-1.5 -1 -0.5 
[Fe/H] 



0.5 



-0.5 



0.5 



-0.5 




1 

0.5 

H 
^ -0.5 
-1 
-1.5 



1 

0.5 

^ 

E -0.5 
-1 H 

-1 .5 



-2.5 



1) 



0.5 



-0.5 




-2.5 -2 -1.5 -1 -0.5 
[Si/H] 




-1.5 -1 -0.5 
[S/H] 




Fig. 17. Observed and predicted abundance ratios versus metallicity for the DLA at Zabs = 2.066 toward Q2231— 00. 
The thick dashed- dotted curves correspond to the "spiral" model at i? = 2 kpc. The thin solid curves correspond to 
the "dwarf irregular" model with a single burst of SF occurring at tb = 0.5 Gyr and with a star formation efficiency 
I' — 4.2 Gyr^^ and a duration Ai = 0.1 Gyr. For the definition of symbols, see Fig. 1131 



reproducing the measured [N/S] and [N/Si] upper Hmits, 
since they predict almost solar [N/S, Si] ratios at [S,Si/H] 
~ —1.25. Such high N/a values at these low metallicities 
have never been observed i n any DLA and H ii region. This 
sugges ts that the yields of Ivan den Hoeck fc GroenwegeiJ 
may lead to an overestimation of the N production 
by intermediate-mass stars at the lower end of masses and 
by low-mass stars, which had time to release N in this DLA 
galaxy, given the inferred long age. In the case of the N 
p roduction by massive intermediate-ma ss stars, the yields 
of lvan den Hoeck fc GroenwegeiJ 1 1997|) are more reliable. 



6.2.3. Q2231-00, z^bs = 2.066 

In this system, the measured abundance ratios at disposal 
are [S/Zn], [S/Fe], [Si/Fe], [Ni/Fe], and the upper limits on 
[N/Si] and [N/S]. We consider the dust-corrected values, 
since non-negligible dust corrections are required in this 
system with [Zn/Fe] — -f0.45±0.07. We first investigated 
the "spiral" model to reproduce the [S/Zn] versus [Zn/H], 
[S/Fe] versus [Fe/H], [Si/Fe] versus [Fe/H] and [Ni/Fe] ver- 



sus [Fe/H] diagrams. The favored "spiral" model turned 
out to be the one representing the inner regions of the 
galactic disk, i.e. at galactocentric radii R < 8 kpc. In the 
inner regions of the disk compared to the outer regions, 
the infall of external gas is faster and leads to higher gas 
densities, and hence to a more efficient SF (see Sect . 16.1 .Til . 
The a-elements thus are more enhanced relative to the Fe- 
peak elements at metallicities [Fe/H] < —1.5 dex in the 
inner regions of the disk than in the outer regions. Despite 
this, the "spiral" model of inner regions of the disk does 
not succeed in reproducing all the measured data points 
(see the thick dashed-dotted curves in Fig. I17|l . It mainly 
fails in correctly reproducing the [S/Zn] versus [Zn/H] and 
[S/Fe] versus [Fe/H] data points, and the [N/Si] versus 
[Si/H] and [N/S] versus [S/H] upper limits. 

We then considered the "dwarf irregular" model. We 
assumed a single burst of SF and explored the follow- 
ing values for the three free parameters characterizing the 
burst: z/ = 2.8 - 5.0 Gyr^^, At = 0.05, 0.1, 0.15 Gyr and 
tb = 0.2, 0.5, 0.8 Gyr. The "dwarf irregular" model which 
best reproduces the data points has one burst with a very 
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high star formation efficiency v = 4.2 Gyr~^ and a dura- 
tion At = 0.1 Gyr (see the thin solid curves in Fig. ITTI). 
The burst occurs after tb — 0.5 Gyr of continuous infall of 
pristine gas. This model is the favored one in comparison 
with the "spiral" model of inner regions of the disk. It re- 
produces at 1 cr all the abundance ratios. Only the [N/S] 
upper limit is at odds. 

We investigated the age of the DLA galaxy assuming 
the constrained "dwarf irregular" model. As in the case 
of the DLA toward QOlOO-l-13, the parameter having the 
highest impact in the age determination is the duration of 
the burst. We feel confident in assuming an age between 
0.05 — 0.15 Gyr, corresponding to the range of values ex- 
plored for At. This yields a redshift of formation between 
2.10 — 2.16 given the assumption that the galaxy has un- 
dergone a single burst of SF. 



6.3. Star formation rates 

The star formation is a key parameter in the formation 
and evolution of galaxies. Therefore, the knowledge of the 
star formation r ate (SFR) of galaxies is very important. 
lKennicufr3 l)l983|) provided the first precise diagnostics for 
the measure of the SFR, such as emission-line fluxes and 
UV contimmm luminosities. Later lMadau et ^ ((l996 h re- 
constructed the cosmic star formation history by measur- 
ing the comoving luminosity density of star-forming galax- 
ies as a func t ion o f redshift. While the original work of 
iMadau et al.l l)l996|) showed a peak in the cosmic star for- 
mation at z ~ 1 — 2, recent results, based on larger sam- 
ples of galaxies and corrections for dust extinction of the 
emitted starlight, do not show such a peak. The SFR per 
unit comoving volume increases by a factor of ~ 10 in the 
redshift interval z = [0,1], and then ei ther remains con- 
stant at z > 1 up to z = 6 llSteidel et al.. ,1999.'l or keens 
increasing to even higher redshifts ((Lanzetta et alJl2002t 
iGalura fc Matteuccil l200,'f^ . However, the galaxies from 
which these results are derived are unlikely to be represen- 
tative of the bulk of the galaxy population in the Universe. 
Indeed, whereas the SFR per unit area for the Milky Way 
is V"^ 4 X 10"^ Mq yr-i kpc"^ (,Kennicutt..l9981 the co- 
moving SFR at ~ 3 is inferred from Lyman-break galax- 
ies, a highly luminous populatio n of star-forming o bjects 
in which ?A > 1 Mq yr"i kpc'^ ||Pettini et al.ll200l|) . As a 
result, the existing measurements of the cosmic star for- 
mation history take into account the contribution of only 
highly luminous and star-forming galaxies. The access to 



the star formation rates in other types of galaxies thus is 
crucial. 

The DLA galaxy population is ideal for this purpose. 
Indeed, these objects sample various types of galaxies over 
a large range of lookback times (see Sects. 16.2. 11 - 1^?^^ . 
since they are detected independently of their distance, 
luminosity and SFH. In addition, because the DLAs are 
not drawn from a flux limited sample of galaxies, we are 
able to derive the SFR values below those determined 
from radiati o n emitted by star-forming galaxies. Recently, 
I Wolfe et al.l 1 2008allJ) provided the first estimations of 
SFRs per unit area in DLAs. Their technique to infer the 
SFRs consists in determining the rate at which the neutral 
gas in DLAs is heated. 

In this work, we also have access to the SFRs of DLAs. 
Indeed, in the case of the "spiral" model the SFR per unit 
area is a direct output of the model, and in the case of the 
"dwarf irregular" model the different physical quantities 
correspond to the absolute values, thus to derive the SFR 
per unit area we assume a spherical symmetry and a galac- 
tic radius of 5 kpc. The derived star formation rates per 
unit area for the DLAs studied are given in Table ITTl and 
are between —2.1 < logip < —1.5 Mq yr""'^ kpc~^. They 
were obtained for the models constrained in the previous 
Sections f Sects. I?).2. 11 - 16. 2.1^ . They correspond in the case 
of the "spiral" model to the SFR that the DLA galaxy has 
at the time of its observations and in the case of the "dwarf 
irregular" model to the average SFR integrated over the 
time of the burst of SF, i.e. over the period when the SF 
is active. The errors on the derived SFRs are very diffi- 
cult to estimate. Our SFR measurements are i n agreement 
with the interval of SFR values obtained bv I Wolfe et alJ 
(2003a b)^ in DLAs using a completely different technique. 
We thus confirm that the SFRs per unit area in DLAs are 
moderate and similar to that measured in the Milky Way 
disk. 

7. Conclusions 

The damped Lya systems are our best laboratory to study 
the high redshift galaxies. Indeed, accurate chemical abun- 
dances of these systems can be obtained over a large in- 
terval of cosmic time, and they thus offer the best oppor- 
tunity to track the chemical evolution of galaxies in the 
Universe. The DLA galaxy population has until now been 



2 The SFRs obtained by IWolfe et al.l ll2003allbl) correspond 
to the SFR values that the DLAs have at the time of their 
observations. 
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analyzed as a whole and chemical evolution models were 
constructed in order to interpret the abundance patterns 
of these objects as an ensemble, considering them as an 
evolutionary sequence. However, several pieces of evidence 
show that the DLAs likely trace galaxies of different types 
and with different evolutionary histories. At this stage of 
knowledge, one would like to define more precisely the star 
formation histories and the chemical evolution stages sam- 
pled by these objects. Therefore, we aimed at constructing 
a sample of DLAs for which one would be able to constrain 
the star formation history, the age and the star formation 
rate of each system individually. For this purpose, it is 
first imperative to obtain comprehensive sets of elemental 
abundances. 

By combining our UVES-VLT spectra of a sample of 
four DLAs in the redshift interval Zabs = 1-7 — 2.5 to- 
ward the quasars Q0100-hl3, Q1331+17, Q2231-00 and 
Q2343-I-12 with the existing HIRES-Keck spectra, we cov- 
ered the total optical spectral range from 3150 to 10 000 
A for the four quasars. Thanks to this large wavelength 
coverage and the high quality of the spectra, we suc- 
ceeded in measuring the column densities of up to 21 
ions, namely of 15 elements — N, O, Mg, Al, Si, P, S, 
CI, Ar, Ti, Cr, Mn, Fe, Ni, Zn. With the detections of 
adjacent ions of the same element, such as A1+/A1++, 
Fe+/Fe++ and N"/N+, and of the Ar/Si,S ratios, we con- 
strained the photoionization effects, which may affect the 
DLA gas-phase abundances. Our analysis revealed that 
the DLA toward Q2343-I-12 requires important ionization 
corrections, while in the three other DLAs the ionization 
corrections are negligible. With the detection of both re- 
fractory and mildly/non-refractory elements of the same 
nucleosynthetic origin, we evaluated the dust depletion cf- 
fe cts and computed the dust corrections using the method 
of lVladild ll2002albl) . The dust corrections are particularly 
important in the DLA toward Q1331-I-17. This system 
with [Zn/Fe] = -|-0.75±0.05 exhibits one of the largest dust 
depletion level of any DLA. The constraint of both the 
photoionization and dust depletion effects allowed us to 
determine relatively robust intrinsic chemical abundance 
patterns of three out of the four DLAs studied. 

The intrinsic chemical abundance patterns are the sig- 
nature of the star formation history of galaxies, which one 
needs to correctly interpret. For this purpose, we called on 
the chemical evo lution models for spiral and dwarf irreg- 
ular galaxies Csee lCalura et al.i.2003'1 . By comparing these 
models with the distributions of the abundance ratios of 
two elements produced on different timescales as a func- 
tion of the metallicity we can constrain the star formation 
history of a DLA galaxy, and with the distributions of the 
same abundance ratios considered as a function of the red- 
shift we can constrain the age of a galaxy. The access to 
several abundance ratios, [S/Zn], [S/Fe], [Si/Fe], [0/Zn], 
[Mg/Fe], [Ni/Fe], [N/S], [N/Si] and [N/0], in the same 
DLA allowed us for the first time to determine the star 
formation history and the age of three DLA systems. We 
were also able to estimate the star formation rates per 
unit area of the DLAs, since the SFR is a direct output of 



the chemical evolution models, once they are constrained. 
Our results show that the galaxy associated with the DLA 
toward QOlOO-l-13 is either a spiral galaxy showing similar 
characteristics as those observed in the solar neighborhood 
(i? ~ 8 kpc) or a dwarf irregular galaxy with a single burst 
of SF with propertiessimilar to the ones observed in the 
BCGs. It is a young galaxy with an age of 250 ± 200 Myr 
or between 50 — 200 Myr, respectively. The galaxy associ- 
ated with the DLA toward Q1331-I-17 is very likely a spiral 
galaxy observed in the outer regions of its disk (i? ~ 16 
kpc) and 1.5 ±0.4 Gyr old. But, it can also well be a dwarf 
irregular galaxy with an inefficient continuous SFH and 
an age > 3.5 Gyr. Finally, the galaxy associated with the 
DLA toward Q2231— 00 is a dwarf irregular galaxy with 
a single intense burst of SF and an age between 50 — 150 
Myr. The estimated star formation rates per unit area of 
these three objects are between — 2.1 < log ijj < — 1.5 Mq 
yr~^ kpc ~^. Thev are in agreement with the values ob- 
tained bv lWolfe et al.l |2003a.h ) in DLAs using a different 
technique, and similar to the SFR per unit area measured 
in the Milky Way disk. 

This work shows that the DLAs studied may either be 
associated with the outer regions of disks of spiral galaxies 
(R > 8 kpc) or with dwarf irregular galaxies with burst- 
ing or continuous star formation histories. Observed at 
redshifts Zabs = 1.7 — 2.5, they may be very young galax- 
ies with ages between 50 — 250 Myr, but also galaxies with 
ages larger than > 1.5 Gyr. Although this work needs to 
be extended to a larger sample of DLA systems to be rep- 
resentative of the overall population of DLAs at high red- 
shift, we can already derive a number of conclusions. Our 
limited sample is centered at z ~ 2, a redshift region where 
spectroscopic data on light-emitting galaxies just start to 
become available. We do not find any association of DLAs 
with episodes of massive star formation, the derived star 
formation rates per unit area are moderate. Clearly this is 
not the standard scenario, and we cannot exclude this pos- 
sibility for all objects in the DLA sample. Actually there 
is an indication of a recent massive burst of star formation 
in at least the DLA at z^hs = 3.025 toward Q0347-3819 
(Lcvshakov et al. 2002; Dcssauges-Zavadskv ct al. 200^). 
The DLAs at z ~ 2 do not appear to be associated with 
a homogeneous population of galaxies, neither in terms of 
star formation history (morphological type) and age. This 
result is in agreement with the studies of DLAs at low red- 
shifts (z < 1) where different types of galaxies have been 
shown to be the optical counterparts (e.g. Brun ct alJ 
llQQTtlNestor et al.ll2002l) . Our results confirm once again 
that simple plots of DLA abundances versus redshift can- 
not be used to estimate the overall metal evolution of the 
Universe, because the local parameters (as determined by 
the star formation history of the associated galaxy) would 
dominate the scatter in the measurements. We plan to ex- 
tend our comparison of the DLA abundance patterns with 
chemical evolution models in 2004 to additional 5 objects 
for which the data have already been acquired. 
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